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Abstract 

We consider solutions to the linear wave equation □gV' = on a suit- 
able globally hyperbolic subset of an extreme Reissner-Nordstrom space- 
time, arising from regular initial data prescribed on a Cauchy hypersurface 
Eo crossing the future event horizon W*". We obtain boundedness, decay, 
non-decay and blow-up results. Our estimates hold up to and including 

. The fundamental new aspect of this problem is the degeneracy of 
the redshift on the event horizon "W^ ■ Several new analytical features of 
degenerate horizons are also presented. 
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D Elliptic Estimates on Lorentzian Manifolds lllll 

1 Introduction 

Black holes are one of the most celebrated predictions of General Relativity, 
and perhaps of Mathematical Physics in general. A particularly interesting but 
peculiar example of a black hole spacetime is given by the so-called extreme 
Reissner- Nordstrom metric, which in local coordinates {t, r, 9, (f) takes the form 

5 = -M2 + ldr2 + r2gg2, (1) 

where 



D = D{r) = 




(7g2 is the standard metric on and M > (see also Appendix [a]). This 
spacetime has been the object of considerable study in the physics literature 
(see for instance the very recent jM!), but the analysis of its waves has not been 
adequately studied from a mathematical point of view. In this paper, we shall 
attempt a more or less complete treatment of the wave equation 

□gV- = (2) 
on extreme Reissner-Nordstrom. Our main results (see Section |4| include: 

1. Local integrated decay of energy, up to and including the event horizon 
%+ (Theorem [l]). 

2. Energy and pointwise uniform boundedness of solutions, up to and includ- 
ing %+ (Theorems [2) [6]). 

3. Energy and pointwise decay of solutions, up to and including (Theo- 
rems [sj [7| § . 

4. Non-decay and blow-up estimates for higher order derivatives of solutions 
along "H"*" (Theorem [9]) . 

The latter blow-up estimates are in sharp contrast with the non-extreme case, 
for which decay holds for all higher order derivatives of tp along H"*". Several 
new analytical features of degenerate event horizons are also presented. The 
above theorems, in particular, resolve Open Problem 4 (for extreme Reissner- 
Nordstrom) from Section 8 of The fundamentally new aspect of this prob- 
lem is the degeneracy of the redshift on the event horizon "H^ . 
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1.1 Preliminaries 



Before we discuss in detail our results, let us present the distinguishing prop- 
erties of extreme Reissner-Nordstrom and put this spacetime in the context of 
previous results. 



1.1.1 Extreme Black Holes 

We briefly describe here the geometry of the horizon of extreme Reissner- 
Nordstrom. (For a nice introduction to the relevant notions, we refer the reader 
to The event horizon corresponds to r = M, where the {t,r) coordi- 

nate system ([T]) breaks down. The coordinate vector field dt, however, extends 
to a regular null Killing vector field T on The integral curves of T on T-l^ 
are in fact affinely parametrized: 

VtT = 0. (3) 

More generally, if an event horizon admits a Killing tangent vector field T for 
which ([3| holds, then the horizon is called degenerate and the black hole extreme. 
In other words, a black hole is called extreme if the surface gravity vanishes on 
H+(see Section 2.3). Under suitable circumstances, the notion of extreme black 



holes can in fact be defined even in case the spacetime does not admit a Killing 
field (see |i7,). 

The extreme Reissner-Nordstrom corresponds to the AI — e subfamily of 
the two parameter Reissner-Nordstrom family with parameters mass M > 
and charge e > 0. It sits between the non-extreme black hole case e < M and 
the so-called naked singularity case M < e. Note that the physical relevance 
of the black hole notion rests in the expectation that black holes are "stable" 
objects in the context of the dynamics of the Cauchy problem for the Einstein 
equations. On the other hand, the so-called weak cosmic censorship conjecture 
suggests that naked singularities are dynamically unstable (see the discussion in 
[23,). That is to say, extreme black holes are expected to have both stable and 
unstable properties; this makes their analysis very interesting and challenging. 



1.1.2 Linear Waves 

The first step in understanding the dynamic stability (or instability) of a space- 
time is by considering the wave equation This is precisely the motivation 
of the present paper. Indeed, to show stability one would have to prove that 
solutions of the wave equation decay sufficiently fast. For potential future appli- 
cations all methods should be robust and the resulting estimates quantitative. 
Robust means that the methods still apply when the background metric is re- 
placed by a nearby metric, and quantitative means that any estimate for ^p must 
be in terms of uniform constants and (weighted Sobolev) norms of the initial 
data. Note also that it is essential to obtain non-degenerate estimates for tp 
on Ji^ and to consider initial data that do not vanish on the horizon. As we 
shall see, the issues at the horizon turn out to be the most challenging part in 
understanding the evolution of waves on extreme Reissner-Nordstrom. 
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1.1.3 Previous Results for Waves on Non- Extreme Black Holes 



The wave equation ^ on black hole spacetimes has long been studied begin- 
ning with the pioneering work of Regge and Wheeler [IT] for Schwarzschild. 
Subsequently, a series of heuristic and numerical arguments were put forth for 
obtaining decay results for ijj (see [SllSn])- However, the first complete quantita- 
tive result (uniform boundedness) was obtained only in 1989 by Kay and Wald 
[52] . extending the restricted result of [5^. Note that the proof of [52] heavily 
depends on the exact symmetries of the Schwarzschild spacetime. 

During the last decade, the wave equation on black hole spacetimes has be- 
come a very active area in mathematical physics. As regards the Schwarzschild 
spacetime, "X estimates" providing local integrated energy decay (see Section 



1.2.1 below) were derived in [6] [71 [18]. Note that [18] introduced a vector field 
estimate which captures in a stable manner the so-called redshift effect, which 
allowed the authors to obtain quantitative pointwise estimates on the horizon 
H+. Refinements for Schwarzschild were achieved in [50] and [3S]. Similar 
estimates to [B] were derived in [5] for the whole parameter range of Reissner- 
Nordstrom including the extreme case. However, these estimates degenerate on 
and require the initial data to be supported away from 

The first boundedness result for solutions of the wave equation on slowly 
rotating Kerr (|a| <C M) spacetimes was proved in [5T] and decay results were 
derived in [S] [551 EH] ■ Decay results for general subextreme Kerr spacetimes 
(|a| < M) are proven in [3S]. Two new methods were presented recently for 
obtaining sharp decay of energy flux and pointwise decay on black hole space- 
times; see [211 [47]. For results on the coupled wave equation see [17]. For other 
results see [221113 IS]- For an exhaustive list of references, see [55] . 

Note that all previous arguments for obtaining boundedness and decay re- 
sults on non-extreme black hole spacetimes near the horizon would break down 



in our case (see Sections 2.3 10.1 13 1. The reason for this is precisely the 



degeneracy of the redshift on 



1.2 Overview of Results and Techniques 

We next describe the results we prove in this paper. All our estimates are with 
respect to norms and we mainly use the robust vector field method. Many 
of our results hold for more general (spherically symmetric) extreme black hole 
spacetimes (but we shall not pursue this here) . We deduce that in some aspects 
the waves on extreme Reissner-Nordstrom have stable properties but in other 
aspects they appear to be unstable; this result is consistent with our discussion 
above about extreme black holes. 



1.2.1 Zeroth Order Morawetz and X Estimates 

Our analysis begins with local spacetime estimates. We refer to local space- 
time estimates controlling the derivatives of ■0 as "X estimates" and ■0 itself 
as "zeroth order Morawetz estimate". Both these types of estimates have a 
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long history (see [55]) beginning with the seminal work of Morawetz [37] for the 
wave equation on Minkowski spacetime. They arise from the spacetime term 
of energy currents associated to a vector field X. For Schwarzschild, such 
estimates appeared in [S] [T] IH] [THl [20] and for Reissner- Nordstrom in [5] . The 
biggest difficulty in deriving an X estimate for black hole spacetimes has to 
do with the trapping effect. Indeed, from a continuity argument one can infer 
the existence of null geodesies which neither cross nor terminate at I"*". In 
our case these geodesies lie on a hypersurface of constant radius (see Section 



2.2 1 known as the photon sphere. From the analytical point of view, trapping 
affects the derivatives tangential to the photon sphere and any non-degenerate 
spacetime estimate must lose (tangential) derivatives (i.e. must require high 
regularity for ip). 

In this paper, we first (making minimal use of the spherical decomposition) 
derive a zeroth order Morawetz estimate for ip which does not degenerate at the 
photon sphere. For the case I > 1 (where I is related to the eigenvalues of the 
spherical Laplacian, see Section [?]) our work is inspired by [20] and for Z = 
we present a method which is robust and uses only geometric properties of the 
domain of outer communications. Our argument applies for a wider class of 
black hole spacetimes and, in particular, it applies for Schwarzschild. Note that 
no unphysical conditions are imposed on the initial data which, in particular, 
are not required to be compactly supported or supported away from Once 
this Morawetz estimate is established, we then show how to derive a degenerate 
(at the photon sphere) X estimate which does not require higher regularity and 
a non-degenerate X estimate (for which we need, however, to commute with the 
Killing T). These estimates, however, degenerate on H"*"; this degeneracy will 
be dropped later. See Theorem [T] of Section |4| 

1.2.2 Uniform Boundedness of Non-Degenerate Energy 

The vector field T = dy is causal and Killing and the energy flux of the current 
jj is non-negative definite (and bounded) but degenerates on the horizon (see 
Section [8|. Moreover, in view of the lack of redshift along H'^ , the divergence of 
the energy current Jf^ associated to the redshift vector field N, first introduced 



in is not positive definite near (see Section 10). For this reason we 



appropriately modify so the new bulk term is non-negative definite near 
y,^. Note that the arising boundary terms can be bounded using Hardy-like 
inequalities. It is important here to mention that a Hardy inequality (in the 
first form presented in Section [g]) allows us to bound the local norm of "0 on 
hypersurfaces crossing "H"*" using the (conserved) degenerate energy of T. Note 
also that the bulk term is not positive far away from and so to control these 
terms we use the X and zeroth order Morawetz estimates. See Theorem [2] of 
Section [l 
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1.2.3 Conservation Laws on ■H+ 

Although the bulk term of the modified redshift current is non-negative definite, 
it degenerates on H'^ with respect to the derivative transversal to This 
degeneracy is a characteristic feature of degenerate event horizons. Indeed, in 



Section 13 we show that the lack of redshift along gives rise to a series of 
conservation laws (see Theorem [s] of Section |4| . These laws apply for waves 
which are supported on the angular frequency I, i.e. for which = for all 
k I, where ipk is the projection of ip (viewed as an function on the spheres 
of symmetry) on the eigenspace 14 of the spherical Laplacian ^ (see Section 
[7|. According to these laws a linear combination of the transversal derivatives 
of i/j of order at most / + 1 is conserved along the null geodesies of H^. As 
we shall see, it is these conserved quantities that allow us to infer the unstable 
properties of extreme black holes, and thus, understanding their structure is 
crucial and essential. A consequence of these laws is that for generic initial 
data, the derivatives transversal to H'^ for waves tp which are supported on 
the low angular frequencies do not decay, and if the order of the derivatives 
is sufficiently high then they in fact blow up along The genericity here 
refers to data for which certain quantities do not vanish on the horizon. 

Based on these conservation laws, we will also prove that the Schwarzschild 
boundedness argument of Kay and Wald [55] cannot be applied in the extreme 
case, i.e. for generic ip, there does not exist a Cauchy hypersurface E crossing 
and a solution ip such that 

TtJj — tJj 

in the causal future of E. 

1.2.4 Sharp Higher Order Estimates 

We next return to the problem of retrieving the derivative tranversal to 
in the X estimate in a neighbourhood of . More generally, we establish 
estimates of the derivatives of ip; see Theorem |4] of Section |4j In view of the 
conservation laws one expects to derive k'th order (fc > 1) estimates close 
to "H^ only if -0; = for all I < k. However, on top of this low frequency 
obstruction comes another new feature of degenerate event horizons. Indeed, to 
obtain such estimates, one needs to require higher regularity for ip and commute 
with the vector field transversal to T-L^ . This shows that exhibits phenom- 



ena characteristic of trapping (see also the discussion in Section 1.3.2). Then 



by using appropriate modifications and the Hardy inequalities we obtain the 



sharpest possible result. See Section 14 Note that although (an appropriate 



modification of) the redshift current can be used as a multiplier for all angular 
frequencies, the redshift vector field N can only be used as a comniutatoi[^ for 
ip supported on the frequencies ^ > 1. 



^The redshift vector field was used as a commutator for the first time in 1211 
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1.2.5 Energy and Pointwise Decay 

Using an adaptation of methods developed in the recent [24] , we obtain energy 
and pointwise decay for tp. See Theorems [5] and [7] of Section [4j In [21], a general 
framework is provided for obtaining decay results. The ingredients necessary 
for applying the framework are: 1) good asymptotics of the metric towards null 
infinity 2) uniform boundedness of energy and 3) integrated local energy decay 
(where the spacetime integral of energy should be controlled by the energy of 
ip and Tip). We first verify that extreme Reissner-Nordstrom satisfies the first 
hypothesis. However, in view of the trapping on the event horizon T^"*", it turns 
out that the method described in [53] can not be directly used to yield decay 
results in the extreme case. Indeed, the third hypothesis of [H] is not satisfied 
in extreme Reissner-Nordstrom. For this reason, we obtain several hierarchies of 
estimates in an appropriate neighbourhood of 'H'^ . These estimates avoid mul- 
tipliers or commutators with weights in t, following the philosophy of [24] . Our 
method applies to black hole spacetimes where trapping is present on (in- 
cluding, in particular, a wide class of extreme black holes). Pointwise estimates 
follow then by commuting with the generators of so(3) and Sobolev estimates. 
We shall also see that the low angular frequencies decay more slowly than the 
higher ones. See Section [TH] 

1.2.6 Higher Order Pointwise Estimates 

In order to provide a complete picture of the behaviour of waves ipj it remains 
to derive pointwise estimates for all derivatives of tp. We show that if ip is 
supported on the angular frequency I, then the derivatives transversal to of 
tp decay if the order is at most I. If the order is / -|- 1, then for generic initial 
data this derivative converges along 7^+ to a non-zero number. As before, by 
generic initial data we mean data for which certain quantities do not vasish on 
H^. If, moreover, the order is at least / -I- 2, then for generic initial data these 
derivatives blow up asymptotically along 'H~^ . On the other hand, for T"^ip one 
needs to consider the order to be at least I + 2 + m so the derivatives blow 
up. See Theorems [s] and [o] of Section |4j Therefore, the T derivativet]^ slightly 
counteract the action of the derivatives transversal to . We conclude this 
paper by deriving similar decay, non-decay and blow-up results for the higher 
order energy. In particular, this shows that one can commute with the redshift 
vector field at most I times for ip supported on the frequency I. See Section [TT] 

1.3 Remarks on the Analysis of Extreme Black Holes 

We conclude this introductory section by discussing several new features of 
degenerate event horizons. 

^It is also shown that Tip decays faster than ip- 
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1.3.1 Dispersion vs Redshift 



In [22] , it was shown that for a wide variety of non-extreme black holes, the red- 
shift on y.'^ suffices to yield uniform boundedness of the non-degenerate energy 
without any need of understanding the dispersion properties of ip. However, in 
the extreme case, the degeneracy of the redshift makes the understanding of the 
dispersion of ip essential even for the problem of boundedness. In particular, 
one has to derive spacetime integral estimates for tp and its derivatives. 



1.3.2 Trapping Effect on 



According to the results of Sections 12 and 14 in order to obtain estimates 
in neighbourhoods of 'H'^ one must require higher regularity for ip and commute 
with the vector field transversal to (which is not Killing). This loss of a 
derivative is characteristic of trapping. Geometrically, this is related to the fact 
that the null generators of H'^ viewed as integrals curves of the Killing vector 
field T are affinely parametrized. 

The trapping properties of the photon sphere have different analytical flavour. 
Indeed, in order to obtain estimates in regions which include the photon 
sphere, one needs to commute with cither T or the generators of the Lie alge- 
bra so(3) (note that all these vector fields are Killing). Only the high angular 
frequencies are trapped on the photon sphere (and for the low frequencies no 
commutation is required) while all the angular frequencies are trapped (in the 
above sense) on 7^+. 



1.3.3 Instability of Degenerate Horizons 

Our setting is appropriate for understading the dynamic formation of extreme 
black holes from gravitational collapse: 



i+ 




S 



One can construct spacetimes (for instance solutions of the Einstein-Maxwell- 
Scalar field system) such that the initial hypersurface 5 ~ K'^ is complete and 
asymptotically flat with one end and the spacetime in region TZ coincides with 
the exterior region of extreme R-N. Consider arbitrary initial data IDsg on Sq 
for the wave equation (in particular, suppose that their support includes H^). 
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By (extending and) solving backwards, we obtain a new initial data set IDs 
on the initial hypersurface S. The initial data /-Dso and IDs are isomorphic. 
Indeed, if one considers a scalar perturbation which corresponds to the data set 
IDs, then the trace of the perturbation on Eg corresponds to the data ID^^. 
Hence, the support of our initial data on Eg should include "H^. That is to say, 
our results represent a poor man's linearisation of the non-linear problem. 
In the linear level, the low angular frequencies not only are they trapped on 
but their evolution is governed by conservation laws along null generators. 
These laws imply that for generic initial data (i.e. data for which certain quan- 
tities do not vanish on H"*" H Eq), the higher order derivatives transversal to 7^+ 
blow up along "H^. These differential operators are translation invariant and do 
not depend on the choice of a coordinate system. The blow-up of these geomet- 
ric quantities suggests that extreme black holes are dynamically unstabl^ We 
hope that the methods we develop in this paper will be useful for proving such 
a result in the nonlinear setting. 

1.4 Open Problems 

An important problem is that of understanding the solutions of the wave equa- 
tion on the extreme Kerr spacetime. This spacetime is not spherically symmetric 
and there is no globally causal Killing field in the domain of outer communi- 
cations (in particular, T becomes spacclike close to the event horizon). Recent 
results (25j overcome these difficulties for the whole subextreme range of Kerr. 
The extreme case remains open. 

One could consider the problem of the wave equation coupled with the 
Einstein-Maxwell equations. Then decay for the scalar field was proven in the 
deep work of Dafermos and Rodnianski 17J. Again these results hold for non- 
extreme black holes. For extreme black holes even boundedness of waves remains 
open. 

1.5 Addendum: Published Version 

After its original appearance on the arxiv, this work was subsequently improved 
with various new results and split into two (now published) parts [21 H] . 

2 Geometry of Extreme Reissner-Nordstrom 
Spacetime 

The extreme Reissner-Nordstrom family is a one parameter subfamily sitting in- 
side the two parameter Reissner-Nordstrom family of 4-dimensional Lorentzian 

•^Price's law conjectures that all "parameters" of the exterior spacetime other than the 
mass, charge and angular momentum — so-caled "hair"- should decay polynomially along the 
event horizon or null infinity. This decay suggests that black holes with regular event horizons 
are stable from the point of view of far away observers and can form dynamicaly in collapse. 
Clearly, our results prove that Price's law does not hold in extreme R-N. 
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manifolds {AfM^e, gM,e), where the parameters are the mass Af > and the 
electromagnetic charge e > 0. The extreme case corresponds to M — e. In 
the ingoing Eddington-Finkelstein coordinates {v, r) the metric of the extreme 
Reissner-Nordstrom takes the form 



g = —Ddv + 2dvdr + r g^2 , 



(4) 



where 



D = D{r)= (I 



M 



and gg2 is the standard metric on The Penrose diagram (see also Appendix 



A. 3) of the spacetime J\f covered by this coordinate system for u e M, r G is 




/ 1- 



We will refer to the hypersurface r = M as the event horizon (and denote 
it by "H^) and the region r < M as the black hole region. The region where 
M < r corresponds to the domain of outer communications. 

In view of the existence of the timelike curvature singularity {r — 0} 'inside' 
the black hole (thought of here as a singular boundary of the black hole region) 
and its unstable behaviour, we are only interested in studying the wave equation 
in the domain of outer communications including the horizon "H^. Note that 
the study of the horizon is of fundamental importance. The horizon determines 
the existence of the black hole and therefore any attempt to prove the nonlinear 
stability of the exterior of black holes must come to terms with the structure of 
the horizon. 

We consider a connected asymptotically flat S0(3)-invariant spacelike hy- 
persurface So in JV terminating at j° with boundary such that SEq = Sq DH^. 
We also assume that if n is its future directed unit normal and T = dy then 
there exist positive constants Ci < C2 such that 

Ci < -.9 ("i n) < C2, 
Ci<-g{n,T) <C2. 

Let M be the domain of dependence of Sq- Then, using the coordinate system 
(v, r) we have 

M = ((-00, +00) X [M , +C50 ) X §2) n J+ (So) , (5) 
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where J+ (Sq) is the causal future of Sq (which by our convention includes 
Eo). Note that 7VI is a manifold with stratified (piecewise smooth) boundary 
dM = (•H+n7W)Ul]o. 




/i- 



2.1 The Foliations and 

We consider the foliation St- = (pri'^o), where ipr is the flow of T = 9„. Of 
course, since T is Killing, the hypersurfaces E7- are all isometric to Sq- 

i+ 




We define the region 

7^(0,r) = Uo<f<rSf. 

On St we have an induced Lie propagated coordinate system (p, w) such that 
p e [A^, +00) and cj e These coordinates are defined such that if Q G E,- 
and Q = {vQ,rQ,ijjQ) then p — tq and w = wq. Our assumption on the normal 
n^o (and thus for ns^) implies that there exists a bounded function gi such 
that 

9p = gidv + 

This defines a coordinate system since [dp^de] — [dp^dfp] — 0. Moreover, the 
volume form of E,- is 

dgT.^ = Vp^dpduj, (6) 
where V is a positive bounded function. 
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One can similarly define an induced coordinate system on S,-. Note that 
only local elliptic estimates (see Appendix [d| are to be applied on St-. 

2.2 The Photon Sphere and Trapping Effect 

One can easily see that there exist orbiting future directed null geodesies, i.e. null 
geodesies that neither cross the horizon 'H+ nor meet null infinity I+. A class 
of such geodesies 7 is of the form 

^->7(r)- (t(r),Q,|,0(T)). 

The conditions V^7 = and g (7, 7) = imply that 

g = 2M, 

which is the radius of the so called photon sphere. The t, 4> depend linearly on 

T. 
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In the subextreme case we have Q = + y 1 ^ ■ fact, from 

any point there is a family of null geodesies which asymptotically converge to 
the photon sphere to the future. The physical interpretation of the photon 
sphere will be crucial in what follows. Indeed, the existence of this "sphere" 
(which is in fact a 3-dimcnsional timelikc hypersurface) implies that the energy 
of some photons is not scattered to null infinity or the black hole region. This is 
the so called trapping effect. Note, in comparison, that in Minkowski spacetime 
all future directed null geodesies meet future null infinity I"*". As we shall see, 
this effect forces us to require higher regularity for the waves in order to achieve 
decay results. 



2.3 The Redshift Effect and Surface Gravity of 

As we have already seen, the vector field dy becomes null on the event horizon 
and is also tangent to it. Therefore, is a null hypersurface. In view of 
the symmetry of the Levi-Civita connection, one easily sees that if = V/, 
where the null hypersurface is given by / = 0, then N generates geodesies and 
satisfies the equation 

Vat AT = AA^. (7) 

A similar equation holds if there exists a Killing field V which is normal to the 
null hupersurface. In this case, we have 

\/yV = kV (8) 

on the hypersurface. Since V is Killing, the function k is constant along the 
integral curves of V. This can be seen by taking the pushforward of Q via the 
flow of V and noting that since the flow of V consists of isometries, the push- 
forward of the Levi-Civita connection is the same connection. The quantity k is 
called the surface gravit]j^oi the null hypersurface. Note that, in Riemannian 
geometry, any Killing field that satisfies ([S]) must have k = 0; this is not the case 
for Lorentzian manifolds, however. It is the equation g (V, V) — that allows 
K not to be zero. In the Reissner-Nordstrom family, the surface gravities of the 

''This plays a significant role in black hole "thermodynamics" (see also [ST] and I40p. 
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two horizons {r = r_} and {r = r^} are given by 



r ± - rzp 
24 



1 dD(r) 

2 dr 



(9) 



where _D is given by ( 86 1 of Appendix |A] Note that in extreme Reissner- 



Nordstrom spacetime we have D (r) = (l — — ) and so r_|_ 



M which 



imphes that the surface gravity vanishes. In general, a horizon whose surface 
gravity vanishes is called degenerate. 

Physically, the surface gravity is related to the so-called redshift effect that 
is observed along (and close to) According to this effect, the wavelength of 
radiation close to H"*" becomes longer as v increases and thus the radiation gets 
less energetic. This effect has a long history in the heuristic analysis of waves 
but only in the last decade has it been used mathematically. For example, 
Price's law (see [T7j) and the stability and instability of Cauchy horizons in 
appropriate setting (see [15]) were proved using heavily this effect. (Note that 
for the latter, one also needs to use the dual blueshift effect which is present 
in the interior of black holes.) Moreover, a second proof of the boundedness 
of waves in Schwarzschild (the first proof was given in [32]) was established in 
[15] based on redshift, whereas in [32] it is proved that under some geometric 
assumptions, the positivity of surface gravity suffices to prove boundedness of 
waves without understanding the trapping. 



3 The Cauchy Problem for the Wave Equation 

We consider solutions of the Cauchy problem of the wave equation ^ with 
initial data 

V-lso = V'o € (So) , ns„Vlso =^16 ^it' (^o) , (10) 

where the hypersurface Sq is as defined in Section [2] and denotes the future 
unit normal of So- In view of the global hyperbolicity of A4, there exists a 
unique solution to the above equation. Moreover, as long as fc > 1, we have 
that for any spacelike hypersurface S 

In this paper we will be interested in the case where k > 2. Moreover, we assume 
that 

lim r^'^ix) = 0. (11) 

For simplicity, from now on, when we say "for all solutions tp of the wave 
equation" we will assume that 1/; satisfies the above conditions. Note 
that for obtaining sharp decay results we will have to consider even higher 
regularity for ip. 
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4 The Main Theorems 



We consider the Cauchy problem for the wave equation (see Section |3| on the 
extreme Reissner-Nordstrom spacetime. This spacetime is partially covered by 



the coordinate systems (f, r), (t, r*), (v, r) and (u, v) described in Appendix A.l 
Recall that M is a positive parameter and D = D{r) = (l — ^) . Recall also 
that the horizon 7^+ is located at {r = M} and the photon sphere at {r — 2M}. 

We denote T = d,j = dt, where dy corresponds to the system (u,r) and dt 
corresponds to (t, r). From now on, dv,dr are the vector fields corresponding 
to {v, r), unless otherwise stated. Note that dr* = dt on and therefore it is 
not transversal to H^, whereas dr is transversal to "H 
The foliations St- and E,- are defined in Sections 



2.1 



and 

associated to the vector field V is defined in Section 5T 
mention that close to 7i~^ we have 
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and the current 



For reference, we 



which degenerates on whereas 

which does not degenerate on T-L^ . As regards the foliation E^., we have for r 
sufficiently large 

where corresponds to the null coordinate system {u,v). The Fourier decom- 
position of on §^ (r) is discussed in Section [Tj where it is also defined what it 
means for a function to be supported on a given range of angular frequencies. 
The notation V'z is also introduced in Section [7] Note that all the integrals are 
considered with respect to the volume form. The initial data are assumed to 
be as in Section |3] and sufficiently regular such that the right hand side of the 
estimates below are all finite. Then we have the following 

Theorem 1. (Morawetz and X estimates) There exists a constant C > 
which depends on M and Eq such that for all solutions ^ of the wave equation 
the following estimates hold 

1. Non-Degenerate Zeroth Order Morawetz Estimate: 

K(0,r) ^ JSo 

2. X Estimate with Degeneracy at and Photon Sphere: 



•R.(0,t) 



r4 




Sn 
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3. X Estimate with Degeneracy at : 

<C f (Jj»).4,+JJ(T*)<J- 

4- Non-Degenerate X Estimate with Commutation: 

If, in addition, tp is supported on the angular frequencies I > 1 then 

Theorem 2. (Uniform boundedness of Non-Degenerate Energy) There 
exists a constant C > which depends on M and Sq such that for all solutions 
tp of the wave equation we have 

Theorem 3. (Conservation Laws along H'^) There exist constants aj,j = 
0, 1, i = 0, 1, .... j + 1, which depend on M and I such that for all solutions 

ip of the wave equation which are supported on the frequency I we have 

j+i 

on H~^. Moreover, there exist constants = 0,1,..., I, which depend on M 
and I such that the quantity 

I 

is conserved along the null geodesies of and thus does not decay for generic 
initial data. 

Theorem 4. (Higher Order Estimates: Trapping on 'H'^) There exists 
ro such that M < Tq < 2M and a constant C > which depends on M , k and 
Eg such that if A = {M < r < ro} then for all solutions ip of the wave equation 
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with tpi = for all I < k — \, fc G N, the following holds 



JU + 
J A 



<cjzj jT^' (TV) n^^^ + cY^f jr° (a;^) 



where Xk = ^ if '4'k = and Xk — otherwise. 

Theorem 5. (Energy Decay) Consider the foliation S,- as defined in Section 
m Let 



and 



where A is as defined in Theorem [7) Here corresponds to the null system 
{u,v). There exists a constant C that depends on the mass M and Eq such that: 

• For all solutions ip of the wave equation we have 



1 



<CE 



where 



• For all solutions ip to the wave equation which are supported on the fre- 
quencies I > 1 we have 

where 
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• For all solutions ip to the wave equation which are supported on the fre- 
quencies I > 2 we have 



where 

E^{iP) = (V^) + (Tt/;) + / {drdrij)n't + / {d,{r^))\ 

Theorem 6. (Pointwise Boundedness) There exists a constant C which 
depends on M and Sq such that for all solutions ip of the wave equation we 
have 

everywhere in TZ, where 

|fe|<2"^^o 

Alternatively, we have 
everywhere in TZ, where 



Theorem 7. (Pointwise Decay) Fix Rq such that M < Rq and let t > 1. 

Let El, E2, E3 be the quantities as defined in Theorem^ Then, there exists a 
constant C that depends on the mass M , Rq and Eq such that: 

• For all solutions ip to the wave equation we have 



\/r ■ t' r ■ ^ 

in {Rq < r}, where 

|fc|<2 

For all solutions ip of the wave equation we have 

r 5 



in {M < r < Rq}, where 

L~(So)- 



Ee ^ Ei+Eiii,)+E4T^)+E5 + Y,{E2{n''i>)+E3{n''^)) + \\drM\'' 



|fe|<2 
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• For all solutions ijj to the wave equation which are supported on the fre- 
quencies I > 1 we have 



3 

r 4 

in {M < r < Rq}, where 

|A:|<2 |fc|<2 

For all solutions tp to the wave equation which are supported on the fre- 
quencies l>2 we have 

T 



in {M <r< Rq}, where 



|fe|<2 



Theorem 8. (Higher Order Estimates I: Decay) Fix Ri such that Ri > M 
and let T > 1. Let also k,l G N. Then there exists a constant C which depend 

on M, /, Ri and Sq such thai the following holds: For all solutions ip of the wave 
equation which are supported on the angular frequencies greater or equal to I, 
there exist norms Ek^uE^^i of the initial data oftp such that 



j;^(5,VK <CEli\forallk<l-2, 



i:^n{M<r<Ri} ' ^ 

• la^Vl < CEki"^ in{M <r< Ri} for all k< 1-2, 

• < CEi_ij\ in{M<r< Ri}, 

T4 

• IS^Vl < CEi^iX in{M<r< 

T4 

Theorem 9. (Higher Order Estimates H: Non-Decay and Blow-up) 

Fix Ri such that Ri > M . Let also fc, Z G N. Then there exists a, positive 
constant c which depends only on M,l,k such that for all solutions ip to the 
wave equation which are supported on the frequency I we have 

• dl+'^tp{T,e,(l)) Hi['il}]{e,(j)) along %+ and generically Hi[tlj]{e,(j)) ^ 
almost everywhere on Sg- 
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\d[+^ij\ (r, e,(j))>c \Hi[i}]{e, 0)1 t'=-i asymptotically onV.+ , for k>2. 
Ju{dri^)n'i — > +00 asT ^ +00 for allk>l + l. 



5 The Vector Field Method 

For understanding the evolution of waves we will use the so-called vector field 
method. This is a geometric and robust method and involves mainly esti- 
mates. The main idea is to construct appropriate (0,1) currents and use Stokes' 
theorem (see Appendix |B]) in appropriate regions. For a nice recent exposition 
see [T]. 

Given a (0,1) current we have the continuity equation 

P,<o = / + I + / V^P^ (12) 

where all the integrals are with respect to the induced volume form and the unit 
normals are future directed. The normal to the horizon "H"*" is also future 
directed and can be chosen arbitrarily. Then the corresponding volume form on 



is defined such that ( 12 ) holds. For definiteness, from now on we consider 
that = dy = T. Furthermore, all the integrals are to be considered with 
respect to the induced volume form and thus we omit writing the measure. 



5.1 The Compatible Currents J, K and the Current 8 

We usually consider currents P^ that depend on the geometry of {A4 , g) and 
are such that both P^ and VP^ depend only on the 1-jet of ijj. This can 
be achieved by using the wave equation to make all second order derivatives 
disappear and end up with something that highly depends on the geometry of 
the spacetime. There is a general method for producing such currents using 
the energy momentum tensor T. Indeed, the Lagrangian structure of the wave 
equation gives us the following energy momentum tensor 

T^. (^) = d^^d^^j - ^g^^d'^^jda^P, (13) 

which is a symmetric divergence free (0,2) tensor. We will in fact consider this 
tensor for general functions : — > M in which case we have the identity 

DivT (V-) = (DgV) d^p- (14) 

Since T is a (0,2) tensor we need to contract it with vector fields of A^. It is 
here where the geometry of Ai makes its appearance. We have the following 
definition 

Definition 5.1. Given a vector field V we define the current by 
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We say that we use the vector field F as a mfiZtjpZier[^ if we apply ( 12 1 for 
the current . The divergence of this current is 

Div( J) = Div (TV) = Div (T) y + T (VV) , (16) 

where (VF)'^' = {g'^'^kVY = . If V is a solution of the wave equation 

then Div (T) = and, therefore, V is an expression of the 1-jet of ip. 

Definition 5.2. Given a vector field V the scalar current is defined by 

K"" (^) = T(v^) (vy) - W {V'vy . (17) 

Note that from the symmetry of the energy momentum tensor T we have 

K"" (^) = r^. (^X^ 

where iTy'' — {Cyg)^'^ is the deformation tensor of V . Clearly if ^/j satisfies the 
wave equation then 

i^^(V.) = v^j;(V'). 

Thus if we use Killing vector fields as multiplier then the divergence vanishes and 
so we obtain a conservation law. This is partly the content of a deep theorem 
of Noethei]^ The first term of the right hand side of ( 16 ) does not vanish when 
we commute the wave equation with a vector field that is not Killing (or, more 
generally, when ^ does satisfy the wave equation). Therefore, we also have the 
following definition. 

Definition 5.3. Given a vector field V we define the scalar current £^ by 

£^ (^) ^DtviT)V^ (Og^) di^ {V) = (DgV) V (V) . (18) 



5.2 The Hyperbolicity of the Wave Equation 

Note that equation ( |13[ ) provides us with an energy momentum tensor for any 
pseudo-Riemannian manifold. However, Lorentzian manifolds admit timelike 
vectors and the hyperbolicity of the wave equation is captured by the following 
proposition 

Proposition 5.2.1. Let Vi, V2 be two future directed timelike vectors. Then the 
quadratic expression T(Vi, V2) is positive definite in d'lj:. By continuity, if one 
of these vectors is null then T(Vi,V2) is non-negative definite in dip. 

Proof. Consider a point p € M and the normal coordinates around this point 
and that without loss of generality V2 = (1, 0, 0, 0). Then the proposition is 
an application of the Cauchy-Schwarz inequality. □ 

An important application of the vector field method and the hyperbolicity 
of the wave equation is the domain of dependence property. 

As we shall see, the exact dependence of T on the derivatives of ip will be 
crucial later. For a general computation see Appendix |A. 2. 2[ 

^the name comes from the fact that the tensor T is muhiphed by V. 

^According to this theorem, any continuous family of isometries gives rise to a conservation 
law. 
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6 Hardy Inequalities 



In this section we establish three Hardy inequahties which as we shah see wiU be 
very crucial for obtaining sharp estimates. They do not require tp to satisfy the 
wave equation and so we only assume that %p satifies the regularity assumptions 
describ ed in Section [S] and (11). From now on, is the foliation introduced in 
Section 



2.1 



(however, these inequalities hold also for the foliation S^). 



Proposition 6.0.2. (First Hardy Inequality) For all functions ip which sat- 
isfy the regularity assumptions of Section we have 



1 

where the constant C depends only on M and Eg 



Proof. Consider the induced coordinate system (p, lo) introduced in Section 2.1 
We use the 1-dimensional identity 

4-C30 /'+00 

{dph)i:'dp = [hi;'] - 2 / hijdpHp 

M JM 

with h — r ~ M . In view of the assumption on t/j, the boundary terms vanish. 
Thus, Cauchy-Schwarz gives 

/ V^dp = -2 / {r - M)^pdptpdp 

J{r>M} J{r>M} 

<2(/ ^P'dp] If {r-Mf{dpi:fdp] . 

\J{r>M} I \J{r>M} I 



Therefore, 



dp< 4 / [r - MY{dpi>Ydp. 

{r>M} J{r>M} 

Integrating this inequality over §^ gives us 

^i^'p'dpdoj < 4 / / D{dptp fp'dpduj. 

§2 J{r>M} P Js^ J{r>M} 



The result follows from the fact that each E^. is diffeomorphic to §^ x [M, 
and the boundedness of the factor V in the volume form ([6]) of E,- . □ 

The importance of the above inequality lies on the weights. The weight of 
{dpipY vanishes to second order on ■H+ but does not degenerate at infinit}!^ 
whereas the weight of ip degenerate at infinity but not at r = M . Similarly, one 



^Note that if we replace r — AI with another function g then we will not be able to make 
this weight degenerate fast enough without obtaining non-trinial boundary terms. 
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might derive estimates for the non-extreme cas^ We also mention that the 
right hand side is bounded by the (conserved) flux of T through (see Section 

§■ 

Proposition 6.0.3. (Second Hardy Inequality) Let tq G {M,2M). Then 
for all functions ijj which satisfy the regularity assumptions of Section^ and any 
positive number e we have 

where the constant depends on M , e, rg and So- 
Proof. We use as before the 1-dimensional identity 

ro _ fro 

M JM 

with h = r — tq. Then 

pro 

(ro - M)i)^{M) = + 2hiljdp%l>dp 



IM 

<e f {dp^Ydp+ I [l + ^]^-'dp, 



ro rro / „2 \ 

2 j„ I / ( 1 I y 1 „;,2, 



M JM 



e 



for any e > 0. By integrating over and noting that the p^ factor that appears 
in the volume form of Sn {r < ro} is bounded we obtain the required result. □ 

The previous two inequalities concern the hypersurfaces Sr that cross H^. 
The next estimate concerns spacetime neighbourhoods of "H^. 

Proposition 6.0.4. (Third Hardy Inequality) Let rQ,ri be such that M < 
To < ri < 2M . We define the regions 

A = 7^(0,T) n {M <r<ro}, 
B = n{0,T) n {ro < r < ri} . 

Then for all functions ip which satisfy the regularity assumptions of Section 
we have 



<c +c D[{d,tpy + {dr^Y 

A Jb JAuB 

where the constant C depends on M , ro, ri and So- 



"These estimates turn out to be stronger since the weight of ip may diverge at r = M in 
an integrable manner. 
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Proof. We again use the l-dimensional identity 



M 



with h such that h = 2{p - M) in [Af, tq] and ft.(ri) = 0. Then 

rro pro rri rri 

2/ V'dp< / / (l-(9p5))^'dp+ / g'(apV)'dp. 

JA/ JA/ Jro JA/ 

Thus, integrating over we obtain 

^^<c f D[{d^i^f + idr^f] + c f ij^, 

for ah f > 0, where the constant C depends on M, tq, ri and Sq. Therefore, 
integrating over f € [0, r] and using coarea formula 



f]df 
'0 / 

completes the proof of the proposition. 



/ 



Uo<T<T 



□ 



7 Elliptic Theory on §^(r),r > 

In view of the symmetries of the spacetime, it is important to understand the 
behaviour of functions on the orbits of the action of S0(3). Clearly, these orbits 
are isometric to S^(r) for r > 0. The most important operator on S^(r) for 
our applications is the spherical Laplacian (i.e. the Laplacian induced by the 
standard metric on S^(r)) which will be denoted by In local coordinates it 
is given by 



1 



deism 684,11;) 



1 



sin 6 



04,34,1!) 



Clearly, /f\ip is a linear second order differential operator 
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For the definition of the Sobolev space see Appendix [c} However, S^(r) is 
compact without boundary and therefore any constant map on the sphere Hes 
in the kernel of /f\ and thus ^ is not invertible. Observing now that 



{-^ + i)i,-^j= / iyv^r+ / (19) 

where y is the induced gradient on the sphere and / the identity map, we 



see that is injective. It is not difficuh to see that (19) imphes that 



O = {—4^ + I) is surjective. Thus we have a well-defined inverse 

0-1 : i^(§2(r)) ^ H\S\r)). 

Using Stokes' theorem on S^(r) it is easy to see that O^^ is self-adjoint. More- 
over, by Rellich theorem, the image i?^(S^(r)) of is compactly imbedded 
in its domain L^(S^(r)). Therefore, is also compact and thus if we view T 
as an endomorphism of L^(S^(r)) then we can apply the spectral theorem and 
so we have a complete orthonormal basis of L^(§^(r)) of eigenvectors of 
with discrete eigenvalues fii \, 0. Hence the numbers — ^ -|- 1 are eigenvalues 
of 4^. By separating the angular variables 0, (j) we can compute explicitly these 
eigenvalues. It turns out that the eigenvalues are equal to e N. The 

dimension of the eigenspaces E'' is equal to 2Z -I- 1 and the corresponding eigen- 
vectors are denoted by y™'',— Z < m < I and called spherical harmonics. We 
have L^(§^(r)) = 'Si>qE'' and, therefore, any function ^p e L^(S^(r)) can be 
written as 



oo / 



"^ = 2^1^ ^m.iY"^''. (20) 

1=0 m=-l 

The right hand side converges to in of the sphere and under stronger 
regularity assumptions the convergence is pointwise. Let us denote by the 
projection of ip onto i.e. 

m— — / 

Proposition 7.0.5. (Poincare inequality) If ^ L'^ (S^(r)) and ipi — for 

all I < L — 1 for some finite natural number L then we have 

and equality holds if and only if tpi =0 for all I ^ L. 
Proof. We have 



1>L 
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Therefore, 



S2(r) 




' 1>L 

where we have used the orthogonaUty of distinct eigenspaces. □ 

This inequaUty is very useful, since it gives us an estimate of the zeroth 
order term. However, it is rather restrictive since it requires tp not to have 
highly symmetric components in its spherical decomposition. This should be 
contrasted with the Hardy inequalities where no such restriction is present. 
However, as we shall see in Section [12] the Poincare inequality is very crucial in 
establishing sharp results. In particular, we will see that the behaviour of the 
spherical decomposition determines the evolution of waves close to H"*" . 

Note also that each eigenspace is finite dimensional (and so complete) 
and so closed. Therefore, L^(S^(r)) — E'- Q) (E'')'^ . For example, we have that 
tjj can be written uniquely as 



>i, 



where ip>i e {E'') = ®i>iE' . 

Returning to our 4-dimensional problem, if ip is sufficiently regular in TZ then 
its restriction at each sphere can be written as in ( |20[ ), where V'm,! = 'Pm,i{v,r) 
and F™^' = Y"^'' {6,(f>). We will not worry about the convergence of the series 
since we may assume that -0 is sufficiently regulai]^ First observe that for each 



summand in ( 20 ) we have 



where S is an operator on the quotient 7W/SO(3). Therefore, if i/j satisfies the 
wave equation then in view of the linear independence of y"^''s the terms 1^171,1 
satisfy S^rn,i = ^^^S^V'm,; and, therefore, each summand also satisfies the wave 
equation. This implies that if ^prn,i = initially then ipm.i = everywhere. 
From now on, we will say that the wave tp is supported on the angular 
frequencies I > L ii t/ji ^ 0,i — 0, ...,L — 1 initially (and thus every- 
where). Similarly, we will also say that t/j is supported on the angular 
frequency I = L ii ip ^ E^ . 

^Indeed, we may work only with functions which are smooth and such that the non-zero 
terms in | |20| are finitely many. Then, since all of our results are quantitative and all the 
constants involved do not depend on ip, by a density argument we may lower the regularity 
of ip requiring only certain norms depending on the initial data of tit to be finite. 
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Another important observation is that the operators dy and dr are endomor- 
phisms of for all /. Indeed, dy commutes with //^ and ii i/; G then 

and thus S^V' S ■ Therefore, if ?A is supported on frequencies I > L then the 
same holds for dri). 



8 The Vector Field T 

One can easily see that — dt — df in the intersection of the corresponding 
coordinate systems. Here dt corresponds to the coordinate basis vector field of 
either {t,r) or (t,r*). Recall that the region M. where we want to understand 
the behaviour of waves is covered by the system {v, r, 9, </>). Therefore we define 
T — dy. It can be easily seen from Q that T is Killing and timelike everywher^^ 
except on the horizon where it is null. 

8.1 Uniform Boundedness of Degenerate Energy 

Recall that 

-fV ± n rp , 

where ttI^ is the deformation tensor of T. Since T is Killing, its deformation 
tensor is zero and so 

= 0. (21) 

Therefore, the divergence identity in the region 7?.(0, r) gives us the following 
conservation law 

J>k + / J>n+ = / (22) 
Since T is null on the horizon H"^ we have that J'^rd^^ > 0. More presicely. 



since n!^_f_ = T, from ([95| (See Appendix A. 2. 2 1 we have 



thus proving the following proposition: 

Proposition 8.1.1. For all solutions ip of the wave equation we have 

Jn^)n'^^ < I Jn^)n^^„- (23) 



^'^Note that in the subextreme range T becomes spacelike in the region bounded by the two 
horizons, which however coincide in the extreme case. 
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We know by Proposition 5.2.1 that Jjn^ is non-negative definite. However, 



we need to know the exact way J^nf^ depends on dip. Note that a;„ (see 
Appendix A. 2. 2 1 is strictly positive and uniformly bounded. However 



which clearly vanishes (to second order) at the horizon. Therefore, ( 95 1 implies 



where the constants in ~ depend on the mass M and Eq. Note that all these 
relations are invariant under the flow of T. 

On H'^, the energy estimate (23) degenerates with respect to the transversal 



derivative drip. It is exactly this that does not allow us to use estimate (23) to 



obtain the boundedness result for the waves in the whole region TZ. However, 
if we restrict our attention to the region where r > tq > M (i.e. away from 



the horizon) then commutating the wave equation with T and estimate ( 23 ) in 
conjunction with elliptic and Sobolev estimates give us the boundedness of ip in 
this region. This result is not satisfactory since it provides no information about 
the behaviour of waves on the horizon (which is the boundary of the black hole) 
and so it is not sufficient for non-linear stability problems. 



9 Integrated Weighted Energy Decay 

It turns out that even proving uniform boundedness of ip requires a strong result 
such as integrated decay. The first result we prove in this direction is that there 
exists a constant C which depends on AI and So such that 

X-(jJ(^X+^2) <C / JJ(V)n^„, 



n 

where the weight x — xi^) degenerates only at H^, at r = 2AI (photon sphere) 
and at infinity. The degeneracy at the photon sphere is expected in view of 
trapping. In particular, as we shall see, such an estimate degenerates only with 
respect to the derivatives tangential to the photon sphere. This degeneracy can 
be overcome at the expense of commuting with T. The degeneracy at is due 
to the lack of redshift and can be overcome by commuting with (the non Killing) 
vector field 9,- and by imposing conditions on the spherical decomposition of ijj 



(see Section 12). In Section 13 we will see that these conditions are necessary. 
The degeneracy at infinity will be dropped in Section [15] where we will use a 
vector field more adapted to the neighbourhoods of (which do not contain 
iO). 

Our approach for obtaining the above estimate is by first deriving the weighted 
estimate for ip and then for the derivatives. Note that no unphysical restric- 
tion on the initial data is required (in particular they are not required to be 
supported away H^). 
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9.1 The Vector Field X 



We are looking for a vector field which gives rise to a current whose divergence is 
non-negative (upon integration on the spheres of symmetry). Our work here is 
inspired by [5D] . We will be working with vector fields of the form X = f {r*) dr- 



(for the coordinate system {t,r*) see Appendix A.l). 



9.1.1 The Spacetime Term 

We compute 

where /' — ^^^jpA and H = | "^^l^^ ■ Note that all the derivatives from now on 
will be considered with respect to r* unless otherwise stated. 

Unfortunately, the trapping obstruction does not allow us to obtain a positive 
definite current K-^ so easily. Indeed, if we assume that all these coefhcients 
are positive then we have 

/'>0,/<0 

fH f f 

- - — > — > -- ^ 
D 2D r 

H 1 

^ r 
H 1 
r 



/>0 
<0, 



but the quantity {jj — changes sign exactly at the radius Q of the pho- 
ton sphere. Therefore, there is no way to make all the above four coefficients 
positive. For simplicity, let us define 

— ^ = -H ■ r + D = . 

9.2 The Case / > 1 

We first consider the case where ^ is supported on the frequencies / > 1. 
9.2.1 The Currents J,f and K^^^ 

To assist in overcoming the obstacle of the photon sphere we shall introduce 
zeroth order terms in order to modify the coefficient of {dtip)"^ and create another 
which is more fiexible. Let us consider the current 

JX^gM,^o = jX ^ g (^) ^y^^ + h (r) {S/ ,w) , 
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where h, w are functions on A4. Then we have 

= /f^ + (V^s) ijV^iP + g (V^VaV) + (V^M V^' (V^w) + 

= /s:^+ .g (V''^Va?A)+ (V^g + 2h\7^'w) V-V^i/; + (V^/iV^w + h [Ugw)) ip^. 

By taking h = l,g = 2G,w — ~G we make the coefficient of t/'V^V' vanish. 
Therefore, let us define 

J^'i = + 2G^ (V^tA) - (V^G) (24) 

Then 

j^x,i ^ ^t^jXA + 2G (VVVaV-) - (OgG) ip^ 

Therefore, if we take G such that 

then 

Note that, since / must be bounded and /' positive, — /"' must become negative 
and therefore has the wrong sign. Note also the factor D at the denominator 
which degenerates at the horizon. One way to overcome this would be to make 
/ approach appropriately the horizon and spatial infinity and /' sufficiently 
concave at the photon sphere. But then, we would need "0^,; = for sufficiently 
large / in order to compensate the loss in the compact intermediate regions. 
However, we want to avoid this restriction on ■0. The way that turns out to 
work is by borrowing from the coefficient of {dr*4')'^ which is accomplished by 
introducing a third current, as in [20] . 

9.2.2 The Current Jff'"^ and Estimates for K^^^ 
We define 
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where j3 will be a function of r* to be defined below. Since 
Dw{fV) = fDw{V) + V{f), 

and 

Div {dr* ) 

we have 



D' 2D 
D r 



D r 



+ 



"4 D 



D'\ f D' _ D" 



D \2r'^ 2D 

Note that the coefficient of ftp^ is independent of the choice of the function (3. 
Let us now take 

„ D X 



r + ' 

where x = r* — a — ^/a and a > a sufficiently large number to be chosen 

appropriately. As wc shall sec, the reason for introducing this shifted coordinate 
X is that we want the origin x = to be far away from the photon sphere. Then 

J_ ^Q^^^ ^ ^^)2 ^ f_P |^^|2 ^ 
1 f" X f" 

4 D 



+ 



a 



a^+x^ D (a2 + x^y D 



D' 



D" 



2r2 2D ■ r 



f 



Wc clearly need to choose a function / that is stricly increasing and changes 
sign at the photon sphere. So, if we choose this function (which we call /") 
such that 

(/")' = ^r^, /" {r* = 0) = /" (r = 2M) = 0, 



then 



P . 1 (/")'' 



X in" 



a 



in' 



1 



a2+x2 D {a^+x^ f D 2D (x^ + a^ f 
If we define X" = f^dr- and 

D' D" 



I = 



2r2 2D ■ r 



r 



then 
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9.2.3 Nonnegativity of K^"''^ 

Proposition 9.2.1. There exists a constant C which depends only on M such 
that for all solutions ip to the wave equation which are supported on the frequen- 
cies l> 1 we have 

I I ^™ ^ < C I K"'-^ »). (26) 



Proof. We have 



K 



For the term / we have 



j^(D' D" \ f^(D^_D^ 
\2r^ 2D-r) ^ 2r \ r D 



However, 



D 



^3 2M AM"^ ^ 
=D2 — - D—, D 



M 
r 

M\ 



I |1 



r / 

J. J y.5 



£>5 





6M 




^3 




6M 




^3 



4M\ 

J.3 j 



Therefore, 



1 = 



r J r° 



Therefore, / > and vanishes (to second order) at the horizon and the photon 
sphere. Note now that the term F is positive x is not in the interval [—a, a]. On 
the other hand, the photon sphere is far away from the region where x € [—a, a] 
so / is positive there. However, / behaves like and so it is not sufficient to 
compensate the negativity of F. That is why we need to borrow from the 
coefficient of J7tjj which behaves like ^ . Indeed, the Poincare inequality gives us 



2 
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and, therefore, it suffices to prove that 



+ F>0 



for all X e [—a, a] or, equivalent ly, r* e [y/a, 2a + y/a\. Then 

„ 1 a;2 — , — 

A, 



3 



D2{x^ + a'^y 2(a;2 + a2) 

with 1 < Aa — >^ 1 as a — >■ +oo. Moreover, since r* > r for big r, by taking a 
sufficiently big we have 

2^—^ > 25. ■' 



{r*f 



with 1 > (5a ^ 1 as a — >■ +oo. Therefore, we need to establish that 

_ (q;2 - a;2) {x + a + y/af Sa 
■~ 4 (x2 + a'^f /« A;^' 

for all a; S [—a, a]. In view of the asymptotic behaviour of the constants A^, 5a 
it suffices to prove that the right hand side of the above inequality is strictly 
less than 1. 

Lemma 9.2.1. Given the function (r*) = (r* (x)) defined above, we have 
the following: If —a<x<Q then 

x + a 
2a2 • 

Also, if Q <x <a then 

p ^ x+a _ ^ 
a;2 + a? 2a 

Proof. For — a < a; < we have 

/" i^) = / l^J^^di > I ^^^^dx > I ^dx '"^ + " 



-V^ + J -a + J -a 2Q:2 20^ 

Now for < a; < a we have 



1 1 \ 

a;^ + a;^ + a;'' + 



1 /""^ /I 1 

dx — I — ^ — -7; ^ 1 



J_Q, x2 + a;2 \a;2 + a2 x'^ + a 

x + a 1 



^ x2 + q2 2o' 



□ 
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Consequently, if —a<x<0 and x — —\a then < A < 1 and 
^ 2(a;2 + a2)3 " 2 (A2 + 1)^ 



3 

2(A2 + 1)^ " ~""3 MO' 



since — > 1 as a — > +cx). Similarly, if < x < a and a; = \a then < A < 1 
and 

a [o? — a;2) (x + a + \Joif a [o? — x^) {x + a + \/a:)'^ 
2(a;2 + a2)^(a2+2ax-a:2) 2(a;2 + a2)^ 

_ (1-A2)(l + Af ^ ^ 
2(1 + A2f ^'^-lO^lO' 

since ^> 1 as a — )• +00. □ 

Note, however, that although the coefficient of ip does not degenerate on the 
photon sphere, the coefficient of the angular derivatives vanishes at the photon 
sphere to second order. Having this estimate for ip, we obtain estimates for its 
derivatives (in Section 9.7 we derive a similar estimate for ip for the case I = 0). 

Proposition 9.2.2. There exists a positive constant C which depends only on 
M such that for all solutions of the wave equation which are supported on the 
frequencies I > 1 we have 



Proof. From ( 26 ) we have that the coefficient of -0 does not degenerate at the 
photon sphere. The weights at infinity are given by the Poincare inequality. 
Commuting the wave equation with T completes the proof of the proposition. 

□ 

9.2.4 The Lagrangian Current 

In order to retrieve the remaining derivatives we consider the Lagrangiarp] 
current 

where 



^^The name Lagrangian comes from the fact that if 1/) satisfies the wave equation and 
C denotes the Lagrangian that corresponds to the wave equation, then C{ip , dip , g~^) = 
gi'-^d^iPd^^iP = Div{5/)V^V) 
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Proposition 9.2.3. There exists a positive constant C which depends only on 
M such that for all solutions ip of the wave equation which are supported on the 
frequencies I > 1 we have 

<j[ fDiv(i;) + cgjr-«--' (TV)V 



Proof. We have 



VD ^,2 VD ,^2 3D fM 



where e > is such that ey D < Therefore, in view of Proposition 
have the required result. 



9.2.2 



we 

□ 



Proposition 9.2.4. There exists a positive constant C which depends only on 
M such that for all solutions -ip of the wave equation which are supported on the 
frequencies I > 1 we have 



Proof. Immediate from Propositions |9.2.1| and |9.2.3[ □ 
9.2.5 The Current Jff"-^ 

In case we allow some degeneracy at the photon sphere we then can obtain 
similar estimates without commuting the wave equation with T. This will be 
very useful whenever our analysis is local. We define the function such that 

in' = ^-^^^ ./v=o)=o. 
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Proposition 9.2.5. There exists a positive constant C which depends only on 
M such that for all solutions ip of the wave equation which are supported on the 
frequencies I > 1 we have 



1 
C 



§2 



where X'^ = f^dr* and the current J^'^ is as defined in Section 
Proof. Recall that 



9.2.1 



K 



x,i 



( fdy j-d p 



D 



AD 



if 



d\lll 



-if 



d\ll 



D' 
D ■ r 



if)' 



( D" 



\2D-r 



iy_ 

2^ 



Note that the coefficient of ijP' vanishes to first order on "H^ (see also Lemma 
9.4. f I and behaves like ^ for large r. Note also that the coefficient of 



converges to M (see again Lemma 9.4.1). The result now follows from Propo- 
sition lOHl □ 

In order to retrieve the 9t-derivative we introduce the current 

L^^ = h'^^V^i^, 

where h'^ is such that: 

For M < r < ro < 2Af, h-" = 

For ro < r < 2M, h'^ < 0, 

For r = 2M, ft,'' = to second order, 

For 2M <r < ri,h'^ < 0, 

For ri < r, h"^ = 

Proposition 9.2.6. There exists a positive constant C which depends only on 
M such that for all solutions tp of the wave equation which are supported on the 
frequencies I > 1 we have 



1 

C 



P • (r - 2M) 



1 



P-ir- 2M) 



< 



§2 



(Div(L;i^) 



r^ 



Proof. We have as before 

Div(L;^') = -^(5tV)' + ^idr'^f + h'' + ^^(^.-v-). 

Since the coefficient of ip{dr*ip) vanishes to first on the horizon, the Cauchy- 
Schwarz inequality and Proposition |9 . 2 . 5| imply the result. □ 
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Finally, we obtain: 



Proposition 9.2.7. There exists a positive constant C which depends only on 
M such that for all solutions of the wave equation which are supported on the 
frequencies I > 1 we have 



1 f fP-{r-2M),^ 1 ,^ ,,2 P ■ {r -2M) ^,2 (r - M) . 



< 



§2 



(Div(<)+i^^^i(V.)+Cif^°'^W). 



Proof. Immediate from Propositions |9.2.1| and |9.2.6[ □ 

9.3 The Case / = 

In the case I = the wave is not trapped. Indeed, if we define 
then we have: 

Proposition 9.3.1. For all spherically symmetric solutions ijj of the wave equa- 
tion we have 

Proof. Immediate from the expression of and the above choice of / = 
f. □ 

9.4 The Boundary Terms 

We now control the boundary terms. 

9.4.1 Estimates for J^ng 

Proposition 9.4.1. Let X — fdr* where f — f{r*) is bounded and S be either 
a 5*0(3) invariant spacelike (that may cross Ti.^ ) or a SO{3) invariant null 
hypersurface. Then there exists a uniform constant C that depends on M , S 
and the function f such that for all we have 



s 
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Proof. Wc work using the coordinate system {v,r,9,ip). First note that 
d dv d dr d 



dr* dr* dv dr* dr 

X^fd„+f-Ddr. 



a„ + Ddr 



Now 



ifn") (d^^py + Dfff [d^ij) {drij) 
+ D 



-Dfn'' - -fn" - -Dfn'' + fvJ 
2 2 2 



{dri^f + 



+ 



-Dfn" ~ -Dfn" - -fn'' 
2 2 2 



The result now foUows from the boundedness of /. □ 
9.4.2 Estimates for J^°^* n^, i = 1, 2 

Proposition 9.4.2. There exists a uniform constant C that depends on M and 
S such that 



where S is as in Proposition \9.4-l\ 
Proof. It suffices to prove 

2G"V'(V^^)-(V;.G")V' + 



(ry 

D 



<B / jJ(V)n^ 



For this we first prove the following lemma that is true only in the case of 
extreme Reissner-Nordstrom (and not in the subextreme range). 

Lemma 9.4.1. The function 



F 



1 



1 



,*\2 



1) 



is bounded inTZU T-L^ . 

Proof. For the tortoise coordinate r* we have 
r*{r) =r + 2M\n{r - M) - 



fr2 



M 

r-M 



C. 
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Clearly, F — ^ as r 

^2 



-oo. Moreover, in a neighbourhood of 'H^ 
1 



(r - My 



Ar < oo. 



This implies the required result. □ 
An immediate collorary of this lemma is that the functions 

Gi = r and 

1 D D 

are bounded in 7?. U V.^. Furthermore, for M <C r we have (/")' ^^,1)^1 
and drD ^ \. Therefore, 



^ r^^ rG"" = — 



2I A/")" rOdrD f^D 



2r 



2r 



2r2 



is bounded. Finally, V^G" = 0. The above bounds and the first Hardy inequal- 
ity complete the proof of the proposition. □ 

9.5 A Degenerate X Estimate 

We first obtain an estimate which does not lose derivatives but degenerates at 
the photon sphere. 

Theorem 9.1. There exists a constant G which depends on AI and Sq such 
that for all solutions ip of the wave equation we have 



n{o,T)\T 



(27) 



Proof. We first decompose V as 

^ = ^>i + V'o 
We apply Stokes' theorem for the current 



in the spacetime region 7?.(0,t) and use Propositions 9.2.7 and 9.4.2 We also 
apply Stokes' theorem in 7?.(0,r) for the current (V'o) and use Proposition 
|9.4.1| and by adding these two estimates we obtain the required result. 

□ 



41 



9.6 A non-Degenerate X Estimate 

We derive an estimate which does not degenerate at the photon sphere but 
requires higher regularity for tp. 

Theorem 9.2. There exists a constant C which depends on AI and So such 
that for all solutions of the wave equation we have 



K(0,- 




(28) 

<c I (jJ(^K^ + jJ(r^KJ. 



Proof. We again decompose ■0 as 



^ = '0>l + IpO 

and apply Stokes' theorem for the current 

J^yi) - (V^>i) + Cjf-^i^>i) + Cjf^\Ti;>,) 



in the spacetime region 7?.(0,r), where C is the constant of Proposition 9.2.4 



and use Propositions 19.2.41 and 19.4.21 Finally, we apply Stokes' theorem in 

,„ , „ , ' — ^ / , s : — ^ i' , I , , ,. 



7?.(0,t) for the current J^' {ipo) and use Proposition 
estimates completes the proof. 



9.4.1 



Adding these two 

□ 



9.7 Zeroth Order Morawetz Estimate for ip 



We now prove weighted estimates of the wave ip itself. In Section 9.2.3 we 
obtained such an estimate for ip>i. Next we derive a similar estimate for the 
zeroth spherical harmonic TpQ. We first prove the following lemma 

Lemma 9.7.1. Fix R > 2M. There exists a constant C which depends on 
M , Sq and R such that for all spherically symmetric solutions ijj of the wave 
equation 



{r=F(.}mi(0,T) 



Proof. Consider the region 

J"(0, r) = {7^(0, r) n {M < r < R}} . 

By applying the vector field X — d^- as a multiplier in the region J^(0, t) we 
obtain 

Jh+ Jt Jt.^ J {r=R}nn{o,T) Jso 



42 




In view of the spatial compactness of the region J- ^ the corresponding space- 



time integral can be estimated using Propositions 9.3.1 and 9.4.1 The boundary 
integrals over Eq and E,- can be estimated using Proposition 9.4. 1[ Moreover, 
for spherically symmetric waves il) we have 



D 



which completes the proof. 
Consider now the region 



g = TZr\{R<r}. 



□ 



Proposition 9.7.1. Fix R > 2M . There exists a constant C which depends on 
M , Eq and R such that for all spherically symmetric solutions ip of the wave 
equation 

1 



r < Ch 



{r=R} 



J; 



Proof. By applying Stokes' theorem for the current J^'^{ip) in the region t/(0, r), 
where again X ~ 9^* (and, therefore, / = 1), we obtain 

JG{0,t) J{r=i?,}ne(0,T) JSf, 

where we have used Proposition |9.4.2| to estimate the boundary integral over 
E^ n g. Note that again rig denotes the unit normal vector to {r = R} pointing 
in the interior of Q. For / = 1 we have 
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where / > and / ^ for r > i? > 2M. If G is the function defined in Section 
then G = Y' ^-^^ therefore for sufficiently large R we have dr*G < 0. 



9.2.1 
TEeS 



-g- 

Since ng — -^dr*, by applying Cauchy- Schwarz for the second term on the 
right hand side (and since the third term is positive) we obtain 



1 

-\ 
e 

for a sufficiently small e. Lemma [9 . 7. 1 1 completes the proof. □ 

It remains to obtain a (weighted) estimate for in the region J-. 

Proposition 9.7.2. Fix sufficiently large R > 2M. Then, there exists a con- 
stant G which depends on M , Eq o,''^d R such that for all spherically symmetric 
solutions ip of the wave equation 

jTi 



T J So 



Proof. By applying Stokes' theorem for the current 
where H = [r — M)^, we obtain 

JJ^(0,t) JdJ^(0,T) 

All the boundary integrals can be estimated using Propositions |9.4.2| and |9.7.1[ 
Note also that 

Since OgH ^ D for M < r < R and ^ is bounded, the result follows in view of 
the spatial compactness of and Proposition |9 . 3 . l] □ 

We finally have the following zeroth order Morawetz estimate: 

Theorem 9.3. There exists a constant C that depends on M and Eq such that 
for all solutions ip of the wave equation we have 

I ^y<C f Jj(^)n^„. (29) 
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Proof. Write 

and use Propositions [9XT| [9J?2] and [Qjl] □ 

Clearly, the lemma used for Proposition |9 . 7. l] holds stricly for the case I = 0. 
In general, one could have argued by averaging the X estimate [27] in R, which 
would imply that there exists a value Rq of r such that all the derivatives are 
controlled on the hypersurface of constant radius Rq. This makes our argument 
work for all ip without recourse to the spherical decomposition. 

Note that the above Morawetz estimate holds for that satisfy the regularity 
assumptions of Section [3] 

Remark 9.1. The key property used in Proposition | . 7.7] is that the d' Alem- 
hertian of ^ is always negative, something not true for larger powers of ^ . Note 
that this unstable behaviour of ^ is expected since it is the static solution of the 
wave equation in Minkowski. 

9.8 Discussion 

In [20], a non-degenerate X estimate is established for Schwarzschild. How- 
ever, there one needs to commute with the generators of the Lie algebra so(3). 
Moreover, the boundary terms could not be controlled by the flux of T but one 
needed a small portion of the redshift estimate. In our case, the structure of r* 
allowed us to bound these terms using only the T flux. 

For a nice exposition of previous work on X estimates see [22j . 

10 The Vector Field N 

It is clear that in order to obtain an estimate for the non-degenerate energy of a 
local observer we need to use timelike multipliers at the horizon. Then uniform 
boundedness of energy would follow provided we can control the spacetime terms 
which arise. For a suitable class of non-degenerate black hole spacetimes, not 
only have the bulk terms the right sign close to but they in fact control the 
non-degenerate energy. Indeed, in [25] the following is proved 

Proposition 10.0.1. Let T-i^ be a Killing horizon with positive surface gravity 
and let V be the Killing vector field tangent to Ti,^ . Then there exists a (f>^ - 
invariant vector field N on and constants b, B > such that for all functions 
tp we have 

onn+. 

The construction of the above vector field does not require the global ex- 
istence of a causal Killing field and the positivity of the surface gravity suf- 
fices. Under suitable circumstances, one can prove that the A'^ flux is uniformly 
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bounded without understanding the structure of the trapping (i.e. no X or 
Morawetz estimate is required). 

However, in our case, in view of the lack of redshift along T-i^ the situation 
is completely different. Indeed, we will see that in extreme Reissner-Nordstrom 
a vector field satisfying the properties of Proposition 1 1 . . T] does not exist. Not 
only will we show that there is no <?!)T--invariant vector field N satisfying the 



properties of Proposition 10.0.1 but we will in fact prove that there is no 



invariant timelike vector field N such that 



K'^ (^) > 



on Ti.'^ . Our resolution to this problem uses appropriate modification of and 
the Hardy inequalities and thus is still robust. 



10.1 The Effect of Vanishing Redshift on Linear Waves 

Let us first try to understand the current K associated to a future directed 
timelike 0T-iiivariant vector field in a neighbourhood of the horizon T-L^ . If 
N = fv (r) dy + fr {r) dr then we obtain 

(V-) =F,v {dy^f + Frr [dr^f + F^ \f^\' + F^r (9.^) (^rV') , 



where the coefficients are given by 



{drfr) _ U 

2 r 



F 



D 



1 



(drfr) , 



Vr 



frD' 



(30) 



Note that since N ~ f^d^ + frdr we have 

g{N,N) = -D{f,f + 2f,fr, 
g{N,T) = -Dfy + fr, 

and so fr {r — M) can not be zero (otherwise the vector field N would not be 



timelike) . Therefore, looking back at the list ( 30 1 we see that the coefficient of 



{dripy^ vanishes on the horizon T-L^ whereas the coefficient of d^^dr^ is equal to 
which is not zero. Therefore, (ip) is linear with respect to dr')p on 
the horizon T-C^ and thus it necessarily fails to be non-negative definite. This lin- 
earity is a characteristic feature of the geometry of the event horizon of extreme 
Reissner-Nordstrom and degenerate black hole spacetimes more generally. This 
proves that a vector field satisfying the properties of Proposition |10.0.1| does 
not exist. 
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10.2 A Locally Non-Negative Spacetime Current 



In view of the above discussion, we need to modify the bulk term by introducing 
new terms that will counteract the presence of dyipdrip. We define 



T _ TN,h jr 



where g is a function on Ai. Then we have 



(32) 



K = K 
= K 
= K 



N,h 



N 



N 



provided ijj is a, solution of the wave equation. Let us suppose that h (r) = "^^^F^ ■ 



Then 



= (^) + h (2d,^dr^ + D [dr^ f + lyVr 

= F,, {d^^f + [Frr + hD] (dr^pf - 



Orfr 



(33) 



Note that by taking h to be constant we managed to have no zeroth order terms 
in the current K. Let us denote the above coefficients of {daipdbip) by Gab where 
a,b € {v, r, y} and define the vector field N in the region M < r < ^ to be 
such that 



and, therefore. 



Mr) = 16r, 
fr{r) = 



M 



(34) 



Clearly, N is timelike future directed vector field. We have the following 

Proposition 10.2.1. For all functions ip , the current K^'^i (ip) defined by (33l 
is non-negative definite in the region A = {M <r< and, in particular, 

there is a positive constant C that depends only on M such that 



K 



Hi^)>C' ((9. V)' + (dr^f + \fM 



Proof. We first observe that the coefficient Gf of dyipdrip is equal to 
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Clearly, G^v = 16 and Grr is non-negative since the factor of the dominant term 
D' (which vanishes to first order on ?^+) is positive. As regards the coefficient 
of the mixed term we have 



Gyr = 16D + 2VD = ei • £2 

where 

62 = 16/D + 2. 
We will show that in region A we have 

£l ^ Grr: 

< Gyy. 

Indeed, we set A = — we have 

' r 
< Grr ^ 

M\ ( M 



1 - 



< 1- 



(1 - A) < (1 - A) 

-A- 



drfrjr) U{r) fr{M) 
2 r M 

3 



+ 



A 



+ A 



which holds. Note also that e\ vanishes to second order on "H^ whereas G 

2 

1 



en < G,,,, 



vanishes to first order. Therefore, in region A we have Grr — 6 

16 

A > 



^ - • 'D. Similarly, 







2 




+ 2 


< 16 ^ 
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which again holds. So in ^ we have 

= a„ + Grr (Pri'f + G™ (9.^9,^) + Gy I^VI^ 

= (G™ - e^) (9,^)2 + (G,, - e?) (9,^)' + 

+ (£2 (5.^))' + (ei + (£2 (ei (9,7^)) + Gy ly^l' 

The above inequalities, the compactness of [Af, ^] and that 

0^+06 + 6^ = ^ (a^ + 5^) + 1 (a + > 
for all a, 6 e E, complete the proof. 



□ 
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10.3 The Cut-off 6 and the Current J^'^' ^ 

Clearly, the current K^'^^ will not be non-negative far away from and thus 
is not useful there. For this reason we extend fy, fr such that 



(r) > for all r>M and U (r) = 1 for all r > 
fr (r) < for all r > M and fr (r) = for all r > 



8M 
8M 



Clearly with the above choices we have g {N, N) < 0,g {N, T) < everywhere 
and so A'' is a future directed timelike (/)^-invariant vector field. 

AT _ 1 . 

Similarly, the modification term in the current ' ^ is not useful for con- 
sideration far away from H^. That is why we introduce a smooth cut-off func- 
tion S : [Af, +oo) R such that S{r) = l,r e [M, ^] and S{r) = 0,r G 
[ : +00 ) . Then we consider the currents 



(35) 



We now consider the three regions A, B, C 











/a / b I 










= and = 





However, this spacetime current, which depends on the 1-jet of '0, will generally 
be negative in region B and thus will be controlled by the X and Morawetz 
estimates (which, of course, are non-degenerate in 6) . 
We next control the Sobolev norm 



IIV'll?ji(E.) 



< 



(9„V)' + (5.Vf + IWl'. 
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Since N and are timelike everywhere in TZ and since uj^ ^uJn are positive 



and uniformly bounded, (95) of Appendix A. 2. 2 implies 

where, in view of the (^y-invariance of Ti^ and N , the constants in ^ depend 
only on M and Sq- Therefore, it suffices to estimate the flux of N through E^. 
For this we have; 

Proposition 10.3.1. There exists a constant C > which depends on M and 
Eg such that for all functions tp 

J^{i^)n^<2[ j;^-'^^hi,)n^ + C I Jl{ij)nr (36) 

S -7- >/ S -7- >J ^Z-Y 



Proof. We have 

for a sufficiently small e. The result follows from the first Hardy inequality. □ 

Corollary 10.1. There exists a constant C > which depends on M and Eq 
such that for all functions ■0 

4''^^h^)n^ <C ( J^{i>)nr (37) 



Proof. Note that 



and use first the Hardy inequality. □ 

10.4 Lower Estimate for an Integral over "H"*" 

Lastly, we need to estimate the integral / J.^^^'- (^)<. ■ Recall that = 
T. We have 
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Proposition 10.4.1. There exist positive constants C and such that for any 
e > and all functions 

Jn+ J'H+ Jt,^ Jt.^ 

where C depends on M , So and depends on M, Eg and e. 
Proof. On we have 5 = 1. Therefore, 

However, 

JH+ J-H+ 

Prom the first and second Hardy inequality we have 
Therefore, 

which completes the proof. □ 

11 Uniform Boundedness of Local Observer's En- 
ergy 

Wc have all tools in place in order to prove the following theorem 

Theorem 11.1. There exists a constant C > which depends on M and Sq 
such that for all solutions ip of the wave equation 

I ^^"'(V'K, < C / J^^(^K,. (38) 

Proof. Stokes' theorem for the current J^''^' ^ in region TZ{0, r) gives us 
>/Et Jn Jh+ Jt.0 
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First observe that the right hand side is controlled by the right hand side of 
(37). As regards the left hand side, the boundary integrals can be estimated 



using Propositions 10.3.1 and 10.4.1 The spacetime term is non-negative (and 
thus has the right sign) in region A , vanishes in region C and can be estimated 
in the spatially compact region B (which does not contain the photon sphere) 
by the X estimate (27) and Morawetz estimate (29). The result follows from 
the boundedness of T-flux through E,- . □ 

Corollary 11.1. There exists a constant C > which depends on M and Eg 
such that for all solutions ip of the wave equation 



K 



(39) 



This corollary and ( 28 1 give us a spacetime integral where the only weight 



that locally degenerates (to first order) is that of the derivative tranversal to 
Recall that in the subextreme Reissner-Nordstrom case there is no such 
degeneration. In the next section, this degeneracy is removed provided ipo = 0- 
This condition is necessary as is shown in Section [T3j 

Corollary 11.2. There exists a constant C > which depends on M and Eq 
such that for all solutions ip of the wave equation 



H + 



(40) 



One application of the above theorem is the following Morawetz estimate 
which does not degenerate at 7^+. 

Proposition 11.0.2. There exists a constant C > which depends on M and 
Eg such that for all solutions ip of the wave equation 



(41) 



Proof. The third Hardy inequality gives us 



AUB 



where C is a uniform positive constant that depends only on M and the regions 
A, B. The integral over B of ip'^ can be estimated using ( 29 ) and the last integral 



on the right hand side can be estimated using (27 1 and Propositions 10.2.1 and 
Corollary [TrT| □ 

Remark 11.1. Note that all the above estimates hold if we replace the foliation 
T,r with the foliation Et- which terminates at 1+ since the only difference is 
a boundary integral over 1+ of the right sign that arises every time we apply 
Stokes ' theorem. The remaining local estimates are exactly th e sa me. However, 
the advantage o/ E^ over E,- will become apparent in Section 15 
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12 Commuting with a Vector Field Transversal 

to n+ 

If we commute the wave equation with T — then we obtain the previous 
estimates for the second order derivatives of tp which involve 9„. Next we 
commute the wave equation □g^/' = with the transversal to the horizon vector 
field dr aiming at controlling all the second derivatives of V' (on the spacelike 
hypersurfaces and the spacetime region up to and including the horizon T-C^). 
Such commutations first appeared in [21j and were used in a more general setting 
for subextreme black holes in |22] . 

12.1 Commutation with the Vector Field dr 

We compute the commutator [Dg , dj] . First note that 

2 
r 

where i? = D' + ^, D' = and 

2 

dr (□g'0) — Ddrdrdrt/j + 2drdydr1p H drdylp + Rdrdrip + drifi^lp 

+ D'drdr^J ~ ^d^tp + R'dri!. 

Therefore, 

□ g [drip) - dr (DgV) = 4^dr1p - V ^ D' drOrtp + \d.,tp - R! dr^p . 



Since 



we obtain 



[/)^,dr]lP=lH. (42) 



[□g, dr] = -D'drdri> + ^d^lp " R' dr^p + ^^^V- 

In case ip solves the wave equation Dg^ = we have 

Ug [drip) = D'drdrip + ^9^,^- - R!drip + ^^V'- (43) 

12.2 The Multiplier L and the Energy Identity 

For any solution -0 of the wave equation we have complete control of the second 
order derivatives of Tp away from "H^ since 

ll'^4'(S.n{M<ro<r<n<2M}) = 1 1 I 1 1 1 (S.n{M<r„ <r}) 

r r (44) 
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where C depends on M and Eq and VV denotes a 2-tensor on the Riemannian 
manifold Y^r and | — | denotes its induced norm (see Appendix |d]). Note that we 
have used the first Hardy inequahty to control the zeroth order terms. Simi- 
larly, away from H^, we have control of the bulk integrals of the second order 
derivatives since 



\dadb 



\L^(TZ(a,T)r]{M<ra<r<ri<2M}) 



(45) 

where C depends on M, rg, ri and Sq and a,b {v, r, J/}. This is proved using 
( 29 ) and local elliptic estimates (appendix [d]) . 

In order to understand the behaviour of the second order derivatives of the 
wave V' in a neighbourhood of T-l~^ we will use the vector field method. We will 
construct an appropriate future directed timelike ^T-invariant vector field 

L = fv9v + frdr, 

which will be used as our multiplier. In view of our discussion above, the vector 
field L will be of interest only in a neighbourhood of the horizon. That is why, L 
will be spatially compactly supported. More precisely, we will construct L such 
that L = in r > ri and L timelike in the region M < r < ri. In particular, we 
will be interested in the region M < r < vq < ri. Note that tq, ri are constants 
which will be determined later on. The regions A, B are depicted below 





i+ 

1 \^ 
i \^ 
1 \^ 


/a / b I 


c 







x+ 



For simplicity we will write 7?. instead of 7?. (0, t), A instead of A (0, r), etc. 
The "energy" identity for the current (5^')/') is 



n 



1^ - 



(46) 

The right hand side is controlled by the initial data and thus bounded. Also 
since L is timelike in the compact region A we have from ( 95 1 of Appendix 
1X2:1 

J^'(9,^AX ^{d^dri^f + {drdr^f + \'fdr^\ (47) 
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where the constants in ~ depend on M, Eo and L. Furthermore, on the horizon 
(where again L is timelike) we have 



(48) 



Therefore, the term that remains to be understood is the bulk integral. Since 
dri^ does not satisfy the wave equation, we have 

= idrij) + ing {dri^))L{drXP). 

We know that 
where 



F.„ 



(drfv) , 

2 r 

\ {drfr) , 

Vr 



D 



I fr 



2 ' 



D {drU 

r 

In view of equation ( |43[ ) we have 

£^ {dri^) - (Dg {dr4'))L{drlP) 

D'drdrip + ^dyi; - R'dri^ + ^^V' 



L (dr^) 



D'fr {drdr^/jf " D' {OM) (9.9.^) ~ R' fv (5.9,-7^) (a,V) 



Jv 



Jr 



2^ (a,9,^) (a,v) + 2^ (9,a,v) (s^v-) - R'U (a.a,^) (^.V') 

+ 2^(d^dr'ip)M + 2^(drdrtp)M. 
r r 

Therefore, we can write 

(dr^b) ^Hi [d^dri^f + H2 [drdr^f + Wdr^]" + 

+ Hi {dydr-ip) (dy^j) + H5 (dydr^j) {dri)) 

+ i/fi (arS^V-) [dv^] + Hj{dydrip)/iiip 

+ H8{drdr'ip)4^ljj + Hq (dydrlp) {drdrlp) + Hio [drdrlp) {drip) 
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where the coefficients Ifj, i = 1, 10 are given by 

(^r/r) fr 



Hi = (drfv) , H2 = D 



r 



3D' 1 

~^fr^ ^3 = - ^ (drfr) , 



Hi = +2^, Hr, = -f^R', He = +2^, Hr = 2^, Hs = 2^, (49) 



H^=D [drU) - D'U - 2y , Hw = -frR'. 

Since L is a future directed timelike vector field we have /„ (r) > and (r) < 
0. By taking drfv [M) sufficiently large we can make Hi positive close to the 
horizon 'H'^ . Also since the term D vanishes on the horizon to second order and 
the terms R, D' to first order and since fr (M) < 0, the coefficient H2 is positive 
close to (and vanishes to first order on it). For the same reason we have 
H^D < and (HgRf < ^ close to %+ . Moreover, by taking ~drfr{M) 
sufficiently large we can also make the coefficient H3 positive close to such 
that Hg < j^. Indeed, it suffices to consider fr such that — "^4?^ ^ _drfrU^ 
and then by continuity we have the previous inequality close to "H"*". Therefore, 
we consider M < < 2M such that in region A ~ {M < r < r^} we have 

/„ > 1. a, A. > 1, -fr > 1, i?l > 1, H2 > 0, i/3 > 1, 

i?^ H-? 9 H9 H^i (50) 

Hs < — , HgD < — , (HgR? Hg<—. 

^ 10' ^ - 10' ^ ^ - 10' ^ 10 

Clearly, tq depends only on M and the precise choice for L close to H'^. In 
order to define L globally, we just extend fy and fr such that 

fv>0 for all r < ri and /„ = for all r > ri, 
—fr > for all r < ri and fr = for all r > ri, 

for some ri such that ro < ri < 2M. Again, n depends only on M (and the 
precise choice for L). Clearly, L depends only on M and thus all the functions 
that involve the components of L depend only on M. 



12.3 Estimates of the Spacetime Integrals 

It suffices to estimate the remaining 7 integrals with coefficients ffj's with i = 
4, 10. Note that all these coefficiens do not vanish on the horizon. We will 

prove that each of these integrals can by estimated by the N flux for ip and Tip 
and a small (epsilon) portion of the good terms in K^{drip). First we prove the 
following propositions. 

Proposition 12.3.1. For all solutions ip of the wave equation and any positive 
number e we have 

I [dr^f < ! eHi{dMf + eH2idrdrtpf 

J A J A 

+cj j;^{i^K^^+cJ Ji^mn^^^, 
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where the constant depends on M , Sq and e. 

Proof. By applying the third Hardy inequaUty for the regions A, B we obtain 

<C f {dr^jf + C f D [{d^dri^f + {drdr^Jjf] , 
JB JAUB 

where the constant C depends on M and Sq. Moreover, since CD vanishes to 
second order at "H^, there exists r^ with M < r^ < such that in the region 
{M < r < r J we have CD < eHi and CD < ei/a- Therefore, 



<C / (dri^)^ + 



{rg <r<ri } 



'{Af<r<rJ 
D [{d^dr^f + [drdr^f] . 



The first integral on the right hand side is estimated using ( 27 1 and the last inte- 
gral using the local elliptic estimate (45 ) since T is timelike in region {r^ < r < ri} 
since M < r^. The result follows from the inclusion {M < r < r^} C A and the 
non-negativity of Hi, i = 1,2 in A. □ 

Proposition 12.3.2. For all solutions ip of the wave equation and any positive 
number e we have 



< a 



where the positive constant C^ depends on M , Sq o'f^ 
Proof. Integrating by parts gives us 



H+ 



1p{dr1p). 



Since ip solves the wave equation, it satisfies 



on H+. Therefore, 



H+ 



1 

M 
1 

2M 



■H+ 



■H + 



2M 



«+ns, 



■H + 



All the integrals on the right hand side can be estimated by 



J^Wn^^ using 



the second Hardy inequality and (40). Furthermore, 



■H+r\T. 
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From the first and second Hardy inequality we have 



J-H+ns JT.0 



where C depends on M and Sq. In addition, from the second Hardy inequality 
we have that for any positive number e there exists a constant which depends 
on M , Eg and e such that 



□ 



Proposition 12.3.3. For all solutions ip of the wave equation and any positive 
number e we have 



< 



SnUE, 



where the positive constant depends on M , Sq and e. 
Proof. In view of the wave equation on T-i^ we have 

2 



■H+ J^-' 



■H + 



■H + 



The first integral on the right hand side can be estimated using Proposition 



12.3.2 For the second integral we have 



H+ 



H + 



where, as above, e is any positive number and Ce depends on M, Eg a-nd e. 



□ 



Estimate for / H4 {d^drip) (dt^ip) 
Jn 



For any e > we have 

i?4 {dvdrip) (dyip) 



where the constant Cg depends only on M,Eo and e. 



58 



Estimate for / Hr, {d^drijj) {d, 
Jn 



As above, for any e > 



A 



J A J A 

J A J A 



where m — max^ e ^H^. Then from Proposition 12.3.1 (where e is replaced 
with — ) we obtain 



H5 (dydrip) (drip) 



J A J A 

+ Cj J^{i^)n>^^+cJ J^imn^^^, 



(51) 



where depends on M, Eg a-nd e. 



Estimate for / Hq (drdrip) {dyijj) 
Jn 



Since 



Divdr — 



Stokes' theorem yields 



■R. 



n 



n 



Sr, is. 



^6 {dv-^P) (dr-ip) dr ■ n-H^ 



In view of the boundedness of the non-degenerate energy we have 

Hr (9,V) {dryj) dr ■ n^, - f H7 (d^ij) [drij] dr ■ n^^ <c([ n^^ 
for some uniform positive constant C which depends on M and Eq . The bound 



ary integral over can be estimated using Proposition 12.3.2 
If Q = drH^ + -Hq, then it remains to estimate 



/ Q (dyiP) (dr^P) + f He (dydr^) (9,V) 
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For the first integral we have 



which can be estimated using (39) and Proposition 12.3.1 The last spacetime 
integral is estimated by (51 1 (where is replaced with H^). 



Estimate for / H-j (d^drip) 
Jn 



n 

Stokes' theorem in region TZ gives us 

Jn 



n 



For the boundary integrals we have the estimate 



< e 



Ans 



where C'^ depends on M, Sq f-nd e. As regards the second bulk integral, after 
applying Stokes' theorem on S^(r) we obtain 



where Cg depends on Af , Eo and e. Note that the second integral on the right 
hand side can be bounded by / J^(TV')n^^. Indeed, we commute with T 

and use Proposition 39 (Another way without having to commute with T is by 
solving with respect to i^ip in the wave equation. As we shall see, this will be 
crucial in obtaining higher order estimates without losing derivatives.) 



Estimate for / Hg (drdrip) /f^ip 
Jn 



We have 



n 



n 



n 



Hs {4^ip) (drtp) dr ■ 



^^8 i4^'P) {drip)dr ■ ns^- 



Hsi^ip) {dr'ijj)dr ■n-H+. 
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The integral over is estimated in Proposition 12.3.3 Furthermore, the 
Cauchy-Schwarz inequality implies 

where e is any positive number and depends only on M, Sq ^-i^d For the 
remaining two spacetime integrals we have 

dr {Hsm (dri^) + f -Hs [H) (dr^) 
A JA^ 

Hs i^^dri' ~ [f , dri^]) (dr^) + f ( OrHs + -H^ (^^) (9,^^) 



(42} 



A J A 

Hs {4^drij) + / (drHs) {/^i^) (dr^) 



r 



= / -Hs lya.^M' + / -drHs (w • fdr^) 

J A J A 



The conditions ( 50 ) assure us that 

and thus the first integral above can be estimated. For the second integral we 
apply Cauchy-Schwarz and take 



-drHs m ■ fdr^) <Ce (V-^^ + ^ / W dri^\\ 

A JT.0 J A 

where again e is any positive number and Ce depends on M , Sq and e. 

Estimate for / H<j {dydrip) (drdrip) 
Jn 

In view of the wave equation we have 

2 r 2 2 



Therefore, 



D Ft 

Hq (dydrlp) {drdri)) = - Hg — {drdr^/jf - Hq — (dri^) ((9^9^-0) 
(dylp) (drdr^) - ^ {/^1p) {drdri>) 

r 2 

Note that in A we have 

HgD<^ 
^ - 10 
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and thus the first term in right hand side poses no problem. Similarly, in A we 
have 



and 



El 

10 



According to what we have proved above, the integral 



(dylp) {drdrip) 



can be estimated (provided we replace iJg with —^). Similarly, we have seen 
that the integral 



can be estimated provided we have Hg < which holds by the definition of 
L. 

Estimate for f Hiq {drdrip) (drip) 

We have 



n 



n 



n 



drHiQ + -HiQ 
r 



Hio (drip) dr ■ nso 
If 5 = drH^ + then 



^5(5,^)2 < max I / (9,^) 



and thus this integral is estimated using Proposition |12.3.T| Also 



^^10 [dr-pf dr ■ - / Hio {dr^f dr ■ riE^ 



for a constant C that depends on M and Eq. Finally we need to estimate the 
integral over H^. It turns out that this integral is the most problematic. In the 
next section we will see the reason for that. Since 

2D' ID 2 



we have 



Hio (M) = -fr (M) R' (M) = -fr (M) ^ > 0. 



(52) 
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In order to estimate the integral along H'^ we apply the Poincarc inequality 



^ {dri>) 

H+ ^ 



< 



■H+ 



Therefore, in view of ( 48 1 it suffices to have 

M2 



1(1 + 1) 



2^a + i)^-(^^)^ 2 



-fr [M] 



2l{l + 1) AAP 

1 1 

< - <^ 

+ - 2 

I > 1. 

Note that for Z = 1 we need to use the whole good term over which appears 
in identity ( 46 ) . This is something we can do, since we have not used this term 
in order to estimate other integrals. That was possible via successive use of the 
Hardy inequalities. 

12.4 Estimates for the Second Order Derivatives 

We can now prove the following theorem 

Theorem 12.1. There exists a constant C > which depends on M and Sq 
such that for all solutions ijj of the wave equation which are supported on the 
frequencies I > 1 we have 



H + 



(53) 



A 



Son^ 



where Xi = */V'' supported on I — \ and Xi = 1 if i' is supported on I >2. 

Proof. First note that L ^ in region A. Thus, the above estimate where Sr 
is replaced by S,- n ^ follows from the estimates ( 44 1 and ( 45 ) , the estimates 
that we derived in Section 12. 3| (by taking e sufficiently small) and the energy 
identity (llel). □ 



Note that the right hand side of ( 53 ) is bounded by 
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This means that equivalently we can apply N as multipher but also as com- 
mutator. One can also obtain similar spacetime estimates for spatially compact 
regions which include A by commuting with T and and applying X as a mul- 
tiplier. Note that we need to commute with in view of the photon sphere. 
Note also that no commutation with the generators of the Lie algebra so(3) is 
required. 



We are now in position to drop the degeneracy of ( 39 1 



Proposition 12.4.1. There exists a uniform constant C > which depends on 
M and Sq such that for all solutions ijj of the wave equation which are supported 
on the frequencies I > 1 we have 

I {dAf <c f J^^(^X„ + c I Jl!{mn^^^ + c f J^^(a,,^x^. 



Proof. Immediate from Proposition 1 12. 3. T] and Theorem |1 2. 1[ □ 

13 Conservation Laws on Degenerate Event Hori- 
zons 

We will prove that the lack of redshift gives rise to conservation laws along 
for translation invariant derivatives. These laws affect the low angular 
frequencies and explain why we need to additionally assume that the zeroth 



spherical harmonic vanishes identically in Proposition 12.4.1 We will also show 
that the argument of Kay and Wald (see [52]) could not have been applied in 
our case. 

We use the regular coordinate system (w,r). Our results are not restricted 
to extreme Reissner-Nordstrom but apply to a very general class of (spheri- 
cally symmetric) degenerate black hole spacetimes. We impose no unphysical 
restriction on the initial data (which, of course, should have finite energy). 

Let us first consider spherically symmetric waves. 

Proposition 13.0.2. For all spherically symmetric solutions "tp to the wave 
equation the quantity 

dr^+j^^ (54) 

is conserved along Ti.^ . Therefore, for generic initial data this quantity does not 
decay. 

Proof. Since ij) solves □gV' — ^^'^ since /f^ij: — we have 

dydr'ip -t- ^dyi! = 

and, since is tangential to V.^ , this implies that drip + jyV' remains constant 
along 7^"*" and clearly for generic initial data it is not equal to zero. □ 
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Proposition 13.0.3. For all solutions to the wave equation that are supported 
on I ~ 1 the quantity 

3 1 

drdrip + J^drtp + J^^J (55) 

is conserved along the null geodesies ofH^. 
Proof. Since /^ip = —^ip, the wave equation on gives us 

2 2 

Moreover, 

2 

r 

and thus by restricting this identity on we take 

29,92^ + ^dydr^ - +^^+(-^+ = 0. (57) 



However, R'{M) = and in view of (56) we have 

o 2 2 2 / 2 

2dyd^ip + —dydrip - Jp9yijj + — I 2dydrip + j^dytp 

which means that the quantity 

3 1 

is constant along the integral curves of T on 7^+. □ 

Proposition 13.0.4. There exist constants a^, i = 0, 1, 2, j = 0, 1 which depend 
on M and I such that for all solutions ip of the wave equation which are supported 
on the (fixed) frequency I, where I >2, we have 



tp = a°dyip + a^dydrip, 
drip = a^dyip + aldydrip + a^dyd^ip, 



onn+. 

Proof. Revisit the proof of Proposition |13.0.3| and use the fact that for alH > 2 
we have -'-^ + R'{M) 7^ 0. □ 

The general form of the conservation laws is the following 
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Theorem 13.1. There exist constants af,j = 0,1,...,/ — I, i — 0, + 1, 

which depend on M and I such that for all solutions ^ of the wave equation 
which are supported on the (fixed) frequency I we have 



i=0 



on T-i^ . Moreover, there exist constants /3i,i — 0,1,...,/, which depend on M 
and I such that the quantity 



1=0 



is conserved along the null geodesies of T-L^ . 

Proof. For the first statement, we will use induction on j for fixed /. According 



to Proposition 13.0.4 the result holds for j = 0, 1. We suppose that it holds for 
j=0,l,...,A: — 1 and we will prove that it holds for j — k provided k < I — 1. 
Clearly, 



1=1 

k 



(58) 



We observe that the coefficients of d^'^'^ip and vanish on Since 

— — '^^Jx'^V': the coefficient of d^^p on is equal to 



2j \lj M2 2 M2 M2 ' 

which is non-zero if and only if / 7^ fc. Therefore, for all fc < / — 1 we solve 
with respect to d^ip and use the inductive hypothesis completing the proof of 



the first statement. For the second statement, we apply (58) for k = I. Then, 
according to our previous calculation, the coefficients of d^ip, and d^ip 

vanish on H^. Therefore, we end up with the terms d^d^ip, fc = 0, 1, / + 1 
and d-l.ip,j = 0, 1, / — 1. Thus, from the first part of the theorem, there exist 
numbers f3i,i ~ 0,1, I which depend on AI and / such that 

/ 

on 'H'^ , which implies that the quantity 

/ 

1=0 
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is constant along the integral curves of T on □ 

An important corollary of the above theorem is the following 

Corollary 13.1. For generic initial data there is no Cauchy hypersurface E 
crossing and a solution ip of the wave equation such that 

Tip — ijj 

in the future of S. 

Proof. Suppose that there exists a wave ip such that Tip = ip. Then we can 
decompose 

V' = -00 + V'>1 

as is shown in Section [T] Then 

Tip =rV'o + T-lpyi 

since T is an endomorphism of the eigenspaces of But Tip = ip and thus 
ipQ = {TiP)q = TipQ. Since ipo is a spherically symmetric wave we have dydripo + 
^dyipo — which yields 

dr-ipo + YT^o = 0. 
M 

However, for generic intial data the above quantity is non zero and this completes 
the proof. □ 

Note that even if we restrict the generic initial data such that the evolving 
waves Ip are supported on / > L, for some L > 1, then the above corollary still 
holds. Indeed, it suffices to decompose 

ip ^1pL + 1p>L + l 

and repeat the argument using again the above theorem. 

This proves that we can not use the argument of Kay and Wald (see [5T] 
and [5 2) ) even for proving the uniform boundedness of the waves. Indeed, this 
argument showed that for any solution of the wave equation ip there is another 
wave Ip such that Tijj = ■0 in the future of a Cauchy hypersurface of the domain 
of outer communications of Schwarzschild. 



14 Sharp Higher Order Estimates 

In this section we commute the wave equation with 9^ where fc e N and k > 2 
aiming at controlling all higher derivatives of ip (on the spacelike hypersurfaces 
and the spacetime region up to and including the horizon H"*"). In view of 



Theorem 13.1 the weakest condition on ip would be such that it is supported on 



the frequencies I > k. 
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14.1 The Commutator [Hg, 
We know that 



therefore, 



Furthermore, 



i=l 
k 



i=l 



Let us compute the commutator [/^, d^] . If we denote 



then 



sin^ 



dg {sm0de) + -^—.d^d^ 



i=0 



k\ 1 



E 



i=0 



Therefore, 



and so 
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14.2 Induction on k 

For any solution ^p of the wave equation we control the higher order derivatives 
in the spacetime region away from 



"'So \i^o J 



where C depends on Af, rg, ri and Eq- This can be proved by commuting 
the wave equation with T*,i = 1, ...,k — 1 and applying X as a multiplier and 
using local elliptic estimates. Thus, we need to understand the behaviour of the 
k^^-oidei derivatives of t/" in a neighbourhood of H^. 

Theorem 14.1. There exists a constant C > which depends on AI , k and Sq 

such that for all solutions ip of the wave equation which are supported on the 
frequencies I > k we have 

A 

1=0 "'So ^^^JsonA 



'So' 



where Xk — ^ if "0 supported on I = k and Xk — ^ if il> is supported on 
I > k + \ . The condition on the spherical decomposition is sharp. 

Proof. We follow an inductive process. We suppose that we have proved that 
there exists a uniform positive constant C that depends on AI, Sq and m such 
that 



J^rDA 



+ 1 aarv-f+x^iyarvi' 

m „ m „ 

,^o"'So j^iJSon^ 



(62) 



^0 



for all natural numbers m < fc — 1 (and thus we have Xm — !)■ Note that for 
m = the above estimate is the result of Section [lO] and for m = 1 of Section 
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In order to derive the above estimate for m = fc we will construct an ap- 
propriate multiplier L*^ = fv^v + fr^r which is a future directed timelike (j)T- 

L'^ — in r > ri 



12 



invariant vector field. Following the spirit of Section 
and timelike in the region M < ri. Again, we will be interested in the region 
A = {M < r < To < ri}. The control for the higher order derivatives will be 
derived from the energy identity of the current 



n 



(63) 

By assumption, the right hand side is bounded. Also since V' is timelike in the 
compact region A we have 



and on the horizon 

It suffices to estimate the bulk integral. We have 
But 

+ i^,,(9„a^v)(9;+V), 



(64) 
(65) 



where the coefficients are given from ([30|. Moreover, equation (60) gives us 
1=1 ^ ^ 1=1 ^ / 



i=l 

k 



i=l 



1=1 

k 



1=1 

k 



Y.[.)f.- diR (9^'+ V) (a.9^v^) - E • • (M'^v) (a.9,^v^) 



j=l 

k 



i=i i=i 
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In order to obtain the sharp result we need the following lemma 

Lemma 14.2.1. Suppose ^ is a solution to the wave equation which is supported 
on the (fixed) frequency I = k. Then for allO < i < k—1 and any positive number 
e we have 

where depends on M , k and Eg. 



Proof. Using Theorem |1 3. 1| it suffices to estimate the integrals 



where < j < fc. For < j < fc — 1 we have 



H+ 

The two boundary integrals can be estimated using the second Hardy inequality. 
As regards the last integral on the right hand side, the Cauchy-Schwarz and 
Poincare inequality imply 

For j = k, we use that {dvd^^){d^^jj) — ^dy{d^ip)^ and the second Hardy 
inequality. □ 

We are now in position to estimate the bulk integrals. We decompose 

ip = ijjk+ ijj>k+i, 

as in Section [7j This is needed in view of the factor Xk in the statement of the 
theorem. 

Estimate for [ if/ ((9,^"'+ V) (^r ^ V) , ^ > 
Jn 

For i = we have 

Hq = -kf^D' > 



and so this coefficient has the "right" sign in (63). 
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For i > 1 we use integration by parts 



From Proposition 12.3.1 applied for d^ip and the inductive hypothesis all 
the integrals over A and Sr can be estimated. It only remains to estimate the 
integral over "H"*" when i — I. In this case, if we apply the Poincare inequality 
for the tpk component we notice that we need to absorb a good term in the 
divergence identity for L'^. Indeed, 



1 A'P 



-fr (M) 



1 M2 

2l{l + l) 



21 [l 



— A- — 
2 7 A/2 + M2 

k 



k 



2 

_.fliM) 
2 



4^ 



This implies that we cannot use the Poincare inequality on "H^ in order to 
estimate {d^4'k)'^ anymore. That is why we proved Lemma 



14.2.1 



which we will 

use for the following integrals. Clearly, for the component tp^k+i we need a 
fraction of [yS^V'!^ along H'^ and thus we can take small (epsilon) portions of 
this term later on. 

Estimate for f {d^d^-'+^'4j) > 1 

For z > 2 we use Stokes' theorem 

In 

-> Jy.^ 
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From Proposition 12.3.1 and the inductive hypothesis we can estimate all the 
above integrals. Note that in order to estimate the integrals along for the 
component V'fc we follow the same argument as in the proof of Lemma |14.2.1[ 
For i = 1 we use the wave equation and thus 



Qk-l 

k-1 



r 



k-1 

3 

k-l 

3 



HldiR (d^-'ij) (a^+V) 



-E 

fe-i 

-E 

The integrals of the first sum can be estimated for j = 0, 1 since their coefficients 
vanish on the horizon and the case j > 2 was investigated above. 

The integrals of the second sum were also estimated before. 

For J = the integral of the third sum can be estimated since its coefficient 
vanishes on 7^+. If j > 1 then again these integrals have been estimated. It 
remains to estimate the integral of the last term. Integration by parts gives 



n 



Hi {d';-^H) (srV) + ^ [drHl + ^Hl ] (arif ^) (a,V) 
Hi (a^^^v^) (a,^^) 



n 



Hi (a^i^&v^) (a,V) a, • n^, - / hI (a^^^v) (a,"v) a, • 



H+ 



Hl{dr^/t^^j) (a,V)a,-n„+. 



If we set Q = dr-Hl + ^Ml then 



k-1 



= j g UarV + E 



n 

k-l 
n 



id: 



,fc-i- 



We estimate this integral by applying Stokes' theorem on §^ and Cauchy- 
Schwarz and using the inductive hypothesis. As regards the last bulk integral 
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we have 



Therefore, by applying Stokes' theorem on S'^ we see that this integral can be 
estimated provided we have 

Hl<~\drfHM). (66) 

Note that does not depend on drfr ■ Furthermore, The boundary integral 
over can be estimated as follows: For the component -0^ we have ^^-i/'fe = 
k{k+i)_^^ aiid_thus d^~^/^ipk depends on d^^/jk for < i < k ~ 1 and thus we 



use Lemma 14.2.1 For the component ?/'>fe+i we commute ^ and we use 
Stokes' theorem on §^ and Cauchy-Schwarz as above. 

Estimate for [ Hf (d^^'^^ij) {d^d^i;),i > 1 
Jn 

The case z = 1 was investigated above. 



For i > 2 we use Proposition 12.3.1 and Cauchy-Schwarz. 

Estimate for [ Hf (d^d^-'ip) {d^d^,ip),i > 1 
Jn 

We use Cauchy-Schwarz and the inductive hypothesis. 

Estimate for [ Hf (^S^^) {dyd^i:),i <k-l 



For i = we solve with respect to in the wave equation and then use 
Cauchy-Schwarz and the inductive hypothesis. We proceed by induction on i. 
Assuming that {/^d-i.ip){dyd^ip) is estimated for all < j < i — 1 we will prove 
that the integral can also be estimated. We have 

The first term on the right hand side can be estimated by solving with respect 
to ^ip in the wave equation and using Cauchy-Schwarz. The second term can 
be estimated by our inductive hypothesis. 

Estimate for / Hf {^sd^-'ip) (d^^^tp) ,i>l 
Jn 
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Integration by parts yields 



-H + 

The second bulk integral and the boundary integrals over E are estimated using 
Stokes' identity on and the inductive hypothesis. The last bulk integral 
is estimated by commuting the spherical Laplacian with 9^ and applying again 
Stokes' identity on and Cauchy-Schwarz. Thus for i > 2 we use the inductive 
hypothesis and for i = 1 it suffices to have 

H'l < -Idrf,' (M) . (67) 

As regards the integral over "H^ we have the following: For the component ipk 
we have 4^ipk — — ^^^s-'^V'fe and thus /^d^~^^pk = — (see Section 

[t]) and so we apply Lemma 14.2.1 to estimate it. For the component '0>fe+i we 



use Stokes' theorem on S^, Cauchy-Schwarz and the inductive hypothesis. 

The construction of L'^ is now clear for all fc S N. It suffices to take (M) < 
0, (M) > and -drf^ and drf^ (M) sufficiently large. As far as the global 
construction for the vector field L'' is concerned, we follow the same idea as in 
Section m 

Note that again no commutation with the generators of the Lie algebra so(3) 
is required. 



The condition on the spherical decomposition is sharp (see Theorem 17.31. 



□ 



15 Energy Decay 

In this section we derive the decay for the non-degenerate energy flux of N 
through an appropriate foliation. Clearly, all the estimates in Section |9] are very 
useful, however they degenerate at infinity. It turns out that one can obtain non 
degenerate estimates on regions which connect 'H'^ and I"*" (without containing 
This has to do with the fact that a portion of the initial energy remains in 
a neighbourhood of Such estimates were first derived in the recent [24l along 
with a new robust method for obtaining decay results. Here we establish several 
estimates which will allow us to adapt the methods of [H] in the extreme case. 
These new estimates are closely related with the trapping properties of H'^. 

We fix i? > 2M and consider the hypersurface Sq which is spacelikc for 
M < r < R and crosses "H^ and for r > i? is given by m = u{po), where po E Eq 
is such that r{po) = R. 
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1+ 




Consider now St- — (^,-(Eo): where (pr is the flow of T. For arbitrary ti < T2 
we define 

~ ^Te[Ti,T2]^T, 

T>ll = nil n {r > i?} , 

Nr = T.rr\{r> R} , 

a;i = nil n {r - i?} . 




15.1 r- Weighted Energy Estimates it a Neighbourhood of 

X+ 

The main idea is to derive a non-degenerate X estimate and then derive similar 
estimates for its boundary terms. From now on we work with the null (u, v) 
coordinates (see Appendix A.l) unless otherwise stated. 

Proposition 15.1.1. Suppose p < 3. There exists a constant C that depends 
on M and Eq such that if ip satisfies the wave equation and cj) — rip then 

I r^^-^+f .-^(p + 2)(^+/ '-^i-pD-rD')m^ 



' '^Tl '''1 



(68) 
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Proof. We first consider the cut-off function C : [R, +00) [0, 1] such that 

C(r) = for all r G 1/2] , 

C(r) = 1 for all r e[R + l, +00) . 

Let q = p — 2. We consider the vector field 

V = rWy 

which we apply as multiplier acting on the function C0 in the region X'^J . Then 



ax's 



For the direction of see Appendix |B] Note that for r > R + 1 we have 
{((!)) = K^i4>) and £^(C0) = f ^(0). Then, 

= 2r«-i(9„0)(a„0) + 9r«-i(a„0)2 + ^ (-gi? - rD') \f<f>\^ . 

Note that since ip solves the wave equation <j) satisfies 

4 D' 

—dudy(p H (f> ~ 4^(f>^0 

D r 

and so 



2 4 

□ g0 = --{du4) - dy4)) - — {dudy(f)) + ^0 

2 £)' 

= --{du<P - 9^,0) H 0, 

r r 

which, as expected, depends only on the 1-jet of </>. Therefore, 

= -2r«-i(9„<^)(c'„<^) + 2r'-i(9,0)2 + D' r'^-^ 4>{d,4>) ■ 

Thus 



However, 



2 



dWd m 



3M 
r 

D'D 

1+^ 
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Note that in Minkowski spacetime we would have no zeroth order term in the 
wave equation. In our case we do have, in such a way, however, such that the 
terms on the right hand side of the above identity have the right sign for p < 3 
and sufficiently largcp^ R. 

In view of the cut-off function C all the integrals over A vanish. Clearly, all 
error terms that arise in the region^ W = supp(C — 1) = {R < r < R+ 1} are 



quadratic forms of the 1-jet of V' and, therefore, in view of (|27|) and (|29|), these 
integrals are bounded by / jj (-0) . Also 



D 



dvll Jn^^ Jn^^ 



4 



The last two integrals on the right hand side appear with the right sign. Finally, 
in view of the first Hardy inequality, the error terms produced by the cut-off C, 
in the region W are controlled by the flux of T through . 

□ 

The reason we introduced the function (j) is because the weight r that it 
contains makes it non-degenerate {ip — Q ox\ 1+ but (f) does not vanish there 



in general). The reason we have divided by in (68) is because we want to 
emphasize the weight that corresponds to ip and not to 

A first application of the above r-weighted energy estimate is the following 

Proposition 15.1.2. There exists a constant C that depends on M and Eq such 
that if ip satisfies the wave equation and as defined above with R sufficiently 
large, then 

(/^ Jj(^)4j dr<cj^ Jj (V) +Cj^ (d^cpf . 

Proof. Applying Proposition |15.1.1| for p — 1 and using the fact that for r > R 
and R large enough 

D-rD' > -, 
2 



we have that there exists a constant C that depends on M and Eq such that 

/ 1 (a„0) V ^ iy</>|' < c / Jjwn| +c [ r-i(a„0)^ (6' 



Note now that since = lyvl , (69) yields 



<C J^Wn^ +C r-'{d,^ 



12 Clearly we need to take R > 2M. 

I'^The weights in r play no role in this region. 



l*This makes this p to be the same as in | 24| . The notation <j> a-^d i/i has however been 
swapped. 
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Furthermore, for sufficiently large R we have 



2 r"^ /t)^2 lU'^r^ 



However, if Q is the cut-off function introduced in the proof of Proposition |15.1.H 
then 

Therefore, the above integral is of the right sign modulo some error terms in the 
region W coming from the cut-off C,. These terms are quadratic in the 1-jet of 



■0 and so can be controlled by estimates ( 27 ) and ( 29 1 . Finally, since is null 



we have >/J(V')"^y ~ (9.0)' + |y0|' and thus by (|69| and the coarea formula 



we have the required result. □ 

This is a spacetime estimate which does not degenerate at infinity. Note 
the importance of the fact that the region does not contain If we are 
to obtain the full decay for the energy, then we need to prove decay for the 
boundary terms in Proposition |15.1.2[ The first step is to derive a spacetime 
estimate of the r-weighted quantity r^^(9.0)^. 

Proposition 15.1.3. There exists a constant C which depends on M and Sq 
such that 



(^J^ r-i(a.0)2^dr <cj^ Jl{^)^t.^+C (9.0)' 



Proof. Appying the r-weighted energy estimate for p = 2 we obtain 



r-^(9.0)2</ ^^|y0|Vc/ J^{i;)nl^ +cf (9.0)' 



The result now follows from ( 69 ) and the coarea formula. □ 



15.2 Integrated Decay of Local (Higher Order) Energy 

We have already seen that in order to obtain a non-degenerate spacetime es- 
timate near T-i^ we need to commute the wave equation with the transversal 
to the horizon vector field dr and assume that the zeroth spherical harmonic 
vanishes. Indeed, if ^ is a spatially compact neighbourhood of "H^ (which may 
contain the photon sphere) then we have: 
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Proposition 15.2.1. There exists a constant C that depends on M and Sq 
such that if -tp satisfies the wave equation and is supported on I > 1, then 

I-T2 



C 



Proof. Immediate from Proposition 12.4.1 (28) and the coarea formula 



□ 



As regards the above boundary terms we have 

Proposition 15.2.2. There exists a constant C that depends on M and Sq 
such that if ^ satisfies the wave equation and is supported on l>2, then 



Proof. Immediate from Theorem 14.1 and the coarea formula. 



□ 



15.3 Weighted Energy Estimates in a Neighbourhood of 

Since for Z = the above estimates do not hold for generic initial data, we are 
left proving decay for the degenerate energy. For this we derive a hierarchy of 
(degenerate) energy estimates in a neighbourhood of . In this section, we 
use the {v,r) coordinates. 

Proposition 15.3.1. There exists a (pr -invariant causal vector field P and a 
constant C which depends only on M such that for all ip we have 

in an appropriate neighbourhood A ofH^. 

Proof. Let our ansatz be P = f^d^ + frdr. Recall that 



where the coefficients are given by 

'{drfr) fA frD' 



D 



1 



{drfr) , 



Fur — ^ (.^rfv) 



Vr 
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Let us take fr{r) = — a/D for AI < r < tq < 2AI, with ro to be determmed 
later. Then 



= D 
= D 

- D. 



D' 
D' 



Vdd' 

2 



4^/D 



the constants in ~ depend on M and the choice for tq. Also, 



K,r = VD 



D{drU) 



< eD 



If we take e sufficiently small and /„ such that 

r,n 2 



D{drU) 



< drfv, 



jT u 



(70) 



then there exists tq > M such that 

in — {M < r < tq}. Extend now P iiiTZ such that fv{r) — 1 and fr{r) — for 
all r > ri > ro for some ri < 2AI. This proves the first part of the proposition. 
In region A we have —g{P,P) > 



D and so ujp ^ VD. Therefore, from (951 



(See Appendix A.2.2| we obtain 



K 



N,5.- 



(71) 



□ 



The vector field P will be very crucial for achieving decay results. It is 
timelike in the domain of outer communications and becomes null on the horizon 
"linearly". This linearity allows P to capture the degenerate redshift in ^ in a 
weaker way than N but in stronger way than T. 

15.4 Decay of Degenerate Energy 
15.4.1 Uniform Boundedness of P-Energy 

First we need to prove that the P-fiux is uniformly bounded. 

Proposition 15.4.1. There exists a constant C that depends on AI and Eq 
such that for all solutions ip of the wave equation we have 



(72) 
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Proof. Stokes' theorem for the current gives us 



1+ J-R 



Note that since P is a future-directed causal vector field, the boundary integrals 
over ^+ and X+ are non-negative. The same goes also for in region A 
whereas it vanishes away from the horizon. In the intermediate region this 



spacetimc integral can be bounded using (27). The result now follows from 



□ 



We are now in a position to derive local integrated decay for the T-cnergy. 

Proposition 15.4.2. There exists a constant C that depends on M and Sq 
such that for all solutions ip of the wave equation we have 



Jj(^)n|^ ]dr<C 



and 



dT<C 



in an appropriate (p^ -invariant neighbourhood A ofH^. 

Proof. From the divergence identity for the current and the boundedness of 
P-energy we have 



<C 

A JSr 



for a uniform constant C. Thus the first estimate follows from (70) and the 



coarea formula. Likewise, the second estimate follows from the divergence iden- 

N S —- 

tity for the current ' ' ^ , the boundedness of the non-degenerate TV-energy 
and (nl. □ 



15.4.2 The Dyadic Sequence p„ 



In view of (28) and Propositions 15.1.2 and 15.4.2 we have 



(73) 



where 
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Moreover, from Propositions 15.4.2 and 15.1.3 we have 



4 Wt < Cli (Tip) + C I J^ {ij)n^ + C I (a„<^)^ (74) 



for a constant C that depends on M and Eg- This imphes that there exists a 
dyadic sequence^ p„ such that 



Pn 



4 W< 



where E'l is equal to the right hand side of ( [74| (with ri = 0) and depends 
only on the initial data of -0. We have now all the tools to derive decay for the 
degenerate energy. 

Proposition 15.4.3. There exists a constant C that depends on M and Sq 
such that for all solutions ip of the wave equation we have 



where Ei is as defined above. 



Proof. We apply (73 1 for the dyadic interval [pn,Pn+i] to obtain 



Pn 



In view of the energy estimate 



JJ(V')< < C Jl{M 



which holds for all t > t' , we have 



(Pn+l - Pn) / JaWn'^ 



< CEi- 



1 



Since there exists a uniform constant b > such that fer„_|_i < r„_|_i — r„ we have 

1 

< CEi- 



-P„+i Pn+l 

Now, for T > pi there exists n e N such that Pn l£ t < Pn+i- Therefore, 



Jj(^)n|^ <C I JJ(V')< 



< 



CEi CE\ 



pI pI+1 



< CEi- 



1 



which is the required decay result for the T-energy. 



□ 



'Dyadic sequence is an increasing sequence p„ such that p„ ^ Pn+i ^ (pn+l — Pn)- 
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15.5 Decay of Non-Degenerate Energy 

We now derive decay for the non-degenerate energy. Note that for obtaining such 
a result we must use Proposition |15.2TT] which however holds for waves supported 
on the frequencies ^ > 1. In this case, in view of the previous estimates we have 

^ ' j;(V^)^^ ]dT <cq (^), (75) 



where 



Proposition 15.5.1. There exists a constant C that depends on Ad and Sq 
such that for all solutions ip to the wave equation which are supported on I > 1 
we have 

where E2 depends only on the initial data of ip and is equal to the right hand 



side of (751 (with ti —0). 



Proof. We apply ( [75| for the interval [0, t] and use the energy estimate 

which holds for all t > t' and for a uniform constant C. □ 
15.5.1 The Dyadic Sequence r„ 



If we consider waves which are supported on Z > 2 then from Propositions 15.1.3 



15.2.2 and by commuting (75) with T we take 



Ex 

Tl 



+ c / j^(a,a,^)n^ +c / (9,0)2 

This implies that there exists a dyadic sequence t„ such that 

li w < ^, 

n I J, 



(76) 



where the constant is equal to the right hand side of (76) (with ti =0). We 
can now derive decay for the non-degenerate energy. 
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Proposition 15.5.2. There exists a constant C that depends on M and Sq 
such that for all solutions ip to the wave equation which are supported on I > 2 
we have 

where is as defined above. 



Proof. If we apply (75 1 for the dyadic intervals [t„,t„ we obtain 

In view of the boundedness of the A^-energy we have 

(r„+i-r„)/ j;(V.K < — ■ 

Since there exists a uniform constant 6 > such that bT„+i < t„+i — t„ we 
obtain 



tN/ i\ ^ CE3 



Now, for T > Ti there exists n G N such that r„ < r < t„-(_i. Therefore, 



./E-r 'n ' n+1 



1+1 

which is the required decay result for the energy. 

Summarizing, we have finally proved the following theorem 



□ 



Theorem 15.1. There exists a constant C that depends on the mass M and 
Eq such that: 

• For all solutions ijj of the wave equation we have 



E, 



J^ii.)nl<CE,^ 



where Ei depends only on the initial data and is equal to the right hand 
side of ([74]) (with ti^O). 

• For all solutions ip to the wave equation which are supported on the fre- 
quencies I > 1 we have 

J^{i^)n% <CE,-, 
E, r 

where E2 depends only on the initial data of ip and is equal to the right 
hand side of (75) (with ti = 0). 
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• For all solutions tp to the wave equation which are supported on I > 2 we 
have 



where depends only on the initial data of tp and is equal to the right 



hand side of {76) ( with ti ~ Q) 



16 Pointwise Estimates 

16.1 Uniform Pointwise Boundedness 

We first prove that all the waves "0 remain uniformly bounded in Ai. We work 
with the foliation (or E,-) and the induced coordinate system {p,uj). For 
To > M we have 



+ 00 



[dpi]]) dp 



< 



ro 

+ C30 



{dpi^fp'dp 



{dpipf p^dp 



Therefore, 



ip'^{ro,Uj)duj < — / / (9pV') p^dpduj 
§2 ra Jg2 

^0 Js^n{r>ro} 



(77) 



where C is a constant that depends only on M and Eg- 

Theorem 16.1. There exists a constant C which depends on M and Eq such 
that for all solutions of the wave equation we have 



r 



(78) 



whe 



^4 = E / Jl^i^' 

|fc|<2'^^0 



Proof. From the Sobolev inequality on §^ we have 



\k\<2 



and the theorem follows from (77) and (38). 



□ 
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^>iI1l~(s,) ^ I |||VV'>i|||i2(s^) + ll|VV?/'>i|||i2(s^) + liin |?A>i 



Note that the above also gives us decay in r. Observe also that for the bound- 
edness result, one could avoid commuting with angular derivatives. Indeed, we 
can decompose 

V' = V'o + V'>i 

and prove that each projection is uniformly bounded by a norm that depends 
only on initial data. Indeed, by applying the following Sobolev inequality on 
the hypersurfaces Tit we have 



(79) 

where again VV is a 2-tensor on the Riemannian manifold and | — | is the 
induced norm. Clearly, C depends only on Sq. Note also that we use the Sobolev 
inequality that does not involve the L^-norms of zeroth order terms. We observe 
that the vector field dy — dr is timelike since g {dy — dr,dy — dr) = —D — 2 
and, therefore, by an elliptic estimate (see Appendix |D]) there exists a uniform 
positive constant C which depends on M and Sq such that 

+ C I j;^{dr^P>i)n^^. 



Pointwise bound for is immidiately derived by (77). Therefore, in view of 



( 38 ) and ( 53 1 we obtain the following: 

There exists a constant C which depends on M and Sq such that for all 
solutions '0 of the wave equation we have 

IV'|'<C/ J^{^)n>^,+Cf J^{mn^^^+C[ J^^(a.VK„. (80) 



16.2 Pointwise Decay 
16.2.1 Decay away %+ 

We consider the region {r > -Ro}, where i?o > M . From now on, C will be a 
constant depending only on M, i?o and Sq- 
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Clearly, in this region we have Jf^n'^ ^ •^J'^s ^ depends on Rq. There- 
fore, from (77) we have that for any r > Rq 



C_ 



Applying the above for "0, fiiVi ^ij^ and using the Sobolev inequality for 
yields 

r ■ 

where C depends on Af, Rq and Sq f^nd 

|fc|<2 

Next we improve the decay with respect to r. Observe that for all r > Rq we 
have 

{ri>f{r,uj)duj^ ( {RQij)\Ro,uj)duj + 2 f f '^dp{pi>)p^dpduj 

§2 JS2 Jg2 J^j, p 



However, from the first Hardy inequality (see Section p| we have 



^ 1.2 ^ r< ( jTi " - 1 



r <c 4 (V')n^^ < CE, 



Moreover, if R is the constant defined in Section 15 and recalling that — 
we have 



/ (dp 

JS^n{r>i?o} 



(pW = / {d,{pi;)r + / (5.0)^ 

Jt^n{R>r>Ro} Jn^ 
J So J No 



88 



where for the second inequality we used Propositions 15.1.1 and 15.1.2 Hence 
for T > 1 we have 



§2 



1 



1 



No 



<CEi-, 



since the quantitiy in the square root is dominated by Ei. Therefore, by the 
Sobolev inequahty on §^ we obtain 



< CE5 



1 



16.2.2 Decay near n+ 

We are now investigating the behaviour of ip in the region {M < r < Rq}. We 
first prove the foUowing 

Lemma 16.2.1. There exists a constant C which depends only on M,Rq such 
that for all ri with M < ri < Rq and all solutions "0 of the wave equation we 
have 



^'^{ri,uj)duj < 



C El 



(n - M)2 



Proof. Using (771 we obtain 
C 



< 



C 



ri Jt^n{r>ri} ri Jj:^n{r>ri} E>[p) 



< 



riD{ri) ys,n{r>ri} 



_C Ei_ 



□ 

Lemma 16.2.2. There exists a constant C which depends only on Af, i?o such 
that for all Tq G [M, Rq], a > and solutions tJj of the wave equation, we have 



1 



iP^{ro,Lj)duj < CEi^^ + CyjE- 



-1 



1? ^ ^ V "'S^n{ro<r<ro+T-°} 



Proof. We consider the hypersurface 7a = {r = ro + T °}. 
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y / \ ^ 

M \ ^ 
/ \ ^ 
/ \ ^ 

/ \ s 








c 




<- / ^ 


II 


So 


/ " 







Then by applying Stokes' theorem for the hypersurfaces shown in the figure 
above we obtain 



r(ro,w)< / ^"(ro + r-",a;) + C' / ^{dp^). 

§2 JS2 JS^n{ro<r<ro+r~"} 



For the first term on the right hand side we apply Lemma f6.2.1 (note that 
M < ro + r~"). The lemma now follows from Cauchy-Schwarz applied for the 
second term, the first Hardy inequality and Theorem |f 5.f [ □ 

The case / = 

We first assume that ij} is spherically symmetric. Then we have the pointwise 
bound 



in {M < r < Rq}, where 



\dr1p\\ 



L~(So) 



This can by easily proved by using the method of characteristics and integrating 
along the char acterist ic u = c the wave equation (expressed in null coordinates). 
Hence Lemma 



16.2.2 



for a = I gives 



< CEq — , 

T5 



where Eg = Ei + Eq. Since ^ is spherically symmetric we obtain 

i^^<CEe\. 

r 5 



(81) 



The case I = 1 
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Suppose that "0 is supported on Z = 1. Then from Lemma 
we obtain 



16.2.2 



for a = J 



g2 T2 T 

1 



< CEi — 5" + Cy El y E2 — 5- 
r 2 



Sxn{ro<r<ro+T-°} 
1 



where we have used Theorem 15.5.2 
Sobolev inequahty on §^ we have 



and £'7 = El + E2. Therefore, by the 



^2 < cEj — 

T2 

in {M <r< Ra}, where Er = E|fe|<2 ^7(f^V)- 
The case I > 2 



Suppose that -0 is supported on Z > 2. Then from (77) and Theorem 15.5.2 



we have that there exists a constant C which depends only on M and Rq such 
that 

in {AI < r < Ra}. By Sobolev we finally obtain 



< CEs 



where 

|fe|<2 

Summarizing, we have proved the following 

Theorem 16.2. Fix Rq such that M < Rq and let r > 1. Let Ei,E2, E3 be the 
quantities as defined in Theorem \15.5.S\ Then there exists a constant C that 
depends on the mass M , Rq and Eq such that: 

• For all solutions ip to the wave equation we have 

\^\^<CE5- ^ 



r ■ 



{Rq < r}, whe 



E,= Y, El (f^V). 

|fc|<2 
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Moreover, 



in {i?o < r}. 

• For all spherically symmetric solutions ip of the wave equation we have 

in{M<r< Ro}, where Eq = E^+ Ei{^) + Ei{T^) + ||a^Vllioo(so)- 

• For all solutions to the wave equation which are supported on the fn 
quency I = 1 we have 



'I' < CE^^ 

r2 



in{M<r< Ro}, where £7 = ^^5 + E|fe|<2 ^2(f^'=V')- 

• For all solutions ip to the wave equation which are supported on the fre 

quencies I >2 we have 

\iP\'<CEs^, 

in {M <r< Ro}, where 



|fe|<2 



17 Higher Order Estimates 

We finish this paper by obtaining energy and pointwise results for all the deriva- 
tives of tp. We first derive decay for the local higher order (non-degenerate) 
energy and then pointwise decay, non-decay and blow-up results. We finally use 
a contradiction argument to obtain non-decay and blow-up results for the local 
higher order energy. 

17.1 Energy and Pointwise Decay 

Theorem 17.1. Fix Ri such that Ri > M and let r > 1. Let also k,l G N. 
Then there exists a constant C which depend on M, /, Ri and So such that the 
following holds: For all solutions ip of the wave equation which are supported 
on the angular frequencies greater or equal to I, there exist norms Eh,i of the 
initial data ofip such that 
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J^(a,V)«| <CEli\foraUk<l-2, 



Jt^n{M<r<Ri} ' ^ ~ 

Jt^n(M<r<RT\ " ^ 



ly^n{M<r<Ri} 

Proof. By commuting with T and applying local elliptic estimates and previous 
decay results, the above integrals decay on n {vq < r < Ri} where tq > M. 
So it sufhces to prove the above result for E,- n A, where ^ is a (/3^^-invariant 
neighbourhood of "H^. For we use the spacetime bound given by Theorem [4] 
which implies that there exists a dyadic sequence r„ such that for all k < I — 1 
we have 



JlUd^i'H <CKi-^ (82) 



where 



Then, by Theorem |4] again we have for any r such that <t < t^+i 

k k 

where E depends only on the initial data. Suppose now that k < 1 — 2. We 
apply Theorem |4] for the dyadic intervals [t„,t„_i] and we obtain 

However, the right hand side has been shown to decay like and thus a 
similar argument as above gives us the improved decay for all A: </ — 2. □ 

The next theorem provides pointwise results for the derivatives transversal 
to H+ of i/j. 

Theorem 17.2. Fix Ri such that Ri > M and let r > 1. Let also k, I, ni e N. 
Then, there exist constants C which depend on M, Z, Ri and Eq such that the 
following holds: For all solutions ip of the wave equation which are supported 
on angular frequencies greater or equal to I, there exist norms E^^i of the initial 
data ip such that 

• \drtp\ < CEk,i- in {M < r < Ri} for allk<l~2, 

• |<9'~V| < CEi^i^i^ m {M <r< Ri}, 

T4, 
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• < CEi^i^ in {M <r< R^}. 



Proof. Let Tq such that M < tq < i?i . We consider the cut-ofF 5 : [M, i?i + 1] ^ 
[0,1] such that 6{r) = 1, for r < Ri + \ and (5(r) = 0, for i?i + l/2 <r < Ri + 1. 

i+ 




In view of Theorem |17.1| if A: < / — 2 then both integrals on the right hand side 
decay hke r^^. If A: = / — 1 then the first integral decays like and the second 
like r"^ and if fc = ^ the first integral decays like r^^ and the second is bounded 
(Theorem |4|. Commuting with the angular momentum operators and using the 
Sobolev inequality yield the required pointwise estimates for d^tp ioi k < I. 

□ 

17.2 Non-Decay 

Let now Hilip] be the function on as defined in Theorem [sj Since Hilip] is 
conserved along the null geodesies of whenever ■0 is supported on the angular 
frequency I, we can simply think of Hi[tp] as a function on §q = Eq H H^. We 
then have the following non-decay result. 

Proposition 17.2.1. For all solutions ip supported on the angular frequency I 
we have 

along T-L^ and generically Hi[ip]{6, (p) ^ almost everywhere on Sq. 
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Proof. Since 



on H'^ and since all the terms in the sum on the left hand side decay (see 



Theorem [lL2| ) we take 6*, 0) Hi[ilj]{e , (j)) on It suffices to show 

that generically Hi[ip]{0, <j>) almost everywhere on Sq. We will in fact show 
that for generic solutions ip of the wave equation the function Hi [ip] is a generic 
eigenfunction of order Z of z& on Sq. 

Note that the initial data prescribed on Eq do not a priori determine the 
function Hiltp] on Sg unless I — 0. Indeed, Hilip] involves derivatives of order 
A; < / + 1 which are not tangential to Sq- For this reason we consider another 
Cauchy problem of the wave equation with initial data prescribed on Sq, where 
the hypersurface Eq is as depicted below: 




Note that the hypersurfaces Eq and Eq coincide for r > Ri. Any initial data 
set prescribed on Eq gives rise to a unique initial data set of Eq and vice versa. 
The Sobolev norms of the initial data on Eq and Eq can be compared using 
the pointwise and energy boundedness. Observe now that given initial data on 
Eg the function Hilip] is completely determined on §q, since Hi[ip] involves only 
tangential to Eg derivatives at Sg. Therefore, generic initial data on Eg give rise 
to generic eigenfunctions Hi of order Z of on Sg. Hence, for generic solutions 
tp of the wave equation the functions Hi [ip] do not vanish almost everywhere on 
§0. 

□ 



17.3 Blow-up 

We next show that the above non-decay results imply that higher order deriva- 
tives of generic solutions blow-up along To make our argument clear we 
first consider the spherically symmetric case where Z = 0. 

Proposition 17.3.1. Let fc G N with k > 2. Then there exists a positive 
constant c which depends only on M such that for all spherically symmetric 
solutions ip to the wave equation we have 

\d^i>\>c\HomT''-' 
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asymptotically on T-i^ . 

Proof. We work inductively. Consider the case k = 2. By differentiating the 



wave equation (see for instance (58l) we take 

2 



(83) 



on Note that Td^ip and drip appear with the same sign. If Ho[ip] = then 
there is nothing to prove. Let's suppose that Hq > 0. Then 



We observe 



Therefore, 



Jo Jo Jo T5 



drip > Hq[iP]t ~ CEeT^- > cHq[iP]t 



iP. 



asymptoticahy on H+. By integrating (83) along H'^ we obtain 

d'ri^ir) = d'rm + jrjdrm - ^^.^(r) - ^^(0) + -^iPir) - ^ 



drip 



1 



1 



M 



M 



1 



< dtm + Y7Am - ( Horn + t^v^ ) - Y^^m + - 



1 



1 



< drm + jjdrm - jjHolip] +CEe^^ 

< -cHq[iP]t 



1 



2M2 



cHo[iP]t 



asymptotically on H^. A similar argument works for any k > 2. Indeed, we 



integrate ([58]) for fc > 1 (and / = 0) along and note that Tdr^'^ip and 
appear with the same sign. Therefore, by induction on fc, the integral /J" drip 
dominates asymptotically all the remaining terms which yields the required 
blow-up rates on T-L'^ . Note that the sign of drip depends on k and Hq[iP]. 

□ 

Corollary 17.1. Let k > 2. For generic initial data which give rise to solutions 
ip of the wave equation we have 



Id^iPl 



-oo 



along 'H^ 
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Proof. Decompose i/" = "00 + '0>i thus 



{M,uj)duj>ATr\d^iPo\ (M,w). 

Hence the resuh follows by commuting with fli, the Sobolev inequahty and the 
fact that the right hand side blows up as t — > +00 as Hq [ip] ^ generically. □ 

Let us consider the case of a general angular frequency I. 

Proposition 17.3.2. Let k,l E N with k > 2. Then there exists a positive 
constant c which depends only on M, I, k such that for all solutions ip to the 
wave equation which are supported on the frequency I we have 



W+''i^\{T,e,ct>)>c\Hime,4>)\ 



asymptotically on 



Proof. We first consider fc = 2. If Hi[iIj]{9, (/)) — then there is nothing to prove. 
Suppose that i//[-0](6', 0) > 0. Note 



Jo Jo : 



i=0 



asymptotically on , since the integral on the right hand side is eventually 



dominated by Hi[tP]t in view of Theorem 17.2 If we integrate ( [58| (applied for 
k = l + l) along the null geodesic of whose projection on the sphere is {6, <f)) 
we will obtain 

5^,+V(t,0» <-ci/;[V'](0,0)T, 

since the integral JJ^ d\^^^ eventually dominates all the remaining terms (again 
in view of the previous decay results). The proposition follows inductively by 
integrating ([58]) as in Proposition |17.3.T Recall finally that for generic solutions 
if) we have Hi [ijj] 7^ almost everywhere on §q . 

□ 

The next theorem provides blow-up results for the higher order non-degenerate 
energy. It also shows that our estimates in Section 14 are in fact sharp (regarding 
at least the restriction on the angular frequencies). 

Theorem 17.3. Fix i?i such that Ri > M. Let also k,l E N. Then for generic 
solutions ■0 of the wave equation which are supported on the (fixed) angular 
frequency I we have 

E^n{Af<r<i?,i} 

as T ^ +00 for all k > I + 1. 
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Proof. Consider M < ro < Ri and let S be the cut-off introduced in the proof 
of Theorem [17^ Then, 

k k+1 

\2 



A n A n 



where C depends on M, Ri and Sq- Then, 

fe 



2 ■ 



where nik G N. Note that for the above inequahty we used local elliptic estimates 
(or a more pedestrian way is to use the wave equation and solve with respect 
to d^ip; this is something we can do since D{ro) > 0). Then using (58) we 
can inductively replace the derivatives with the T derivatives. Therefore, 
commuting with Qi and applying the Sobolev inequality imply that for any 
ro > M we have \d^ip\ — as r +oo along r = tq. Let us assume now 
that the energy of d^ip on E,-^ n {M < r < _Ri} is uniformly bounded by B (as 
Tj — !■ +CX)). Given e > take vq such that tq — M — and let be such 

that for all r > we have |(9^-0(r, tq) | < Then, 

|9^-V(r,,M)| < / |9^-V(r„ro)|+ / \dpd^i:\ 

/S2 Jl:-r.n{AI<r<ra} 



2 




<^+ro(ro-Af)M / Jj! (d^i^H <^ 



^.n{M<r<B^i} 



for all T > This proves that J^2 ^)| ~^ *-* ^ ~^ along 'H'^. 

However, in view of Propositions |17.2.1| and |17.3.2| we have 

We have seen that for generic ip the function Hi [ip] is non-zero almost everywhere 
and since it is smooth we have J^2(^j^j^ > 0. This shows that the integral 

/g2 \d^^p{Tj,M)\ can not decay, contradiction. □ 
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A On the Geometry of Reissner-Nordstrom 



The coupled Einstein-Maxwell equations consist of the system 

V^F^, = 0, 
dF = 0, 

where g is Lorentzian metric on an appropriate manifold Ai and (g) , R (g) 
are the Ricci and scalar curvature of the Levi-Civita connection, respectively 
and F a 2- form on Ai . Here T^^ denotes the electromagnetic energy momentum 
tensor. 

The unique family of spherically symmetric asymptotically flat solutions of 
these equations is the two parameter Reissner-Nordstrom family of 4-dimensional 
Lorentzian manifolds {Mm.ct gM.e) where the parameters M and e are called 
mass and (electromagnetic) charge, respectively. The extreme case corresponds 
to M = |e|. 



A.l Constructing the Extention of Reissner-Nordstrom 

We first present the Reissner-Nordstrom metric in local coordinates {t, r, 9, (f) 
which were discovered in 1916 |32] and 1918 [31]. In these coordinates, one 
metric component blows up for various values of r and it is not a priori obvi- 
ous what is the appropriate underlying manifold to study the geometry of this 
solution. Indeed, as we shall see, one can construct another coordinate system 
(u, r, 9, (p) which covers a "bigger" manifold (which we will denote by M.) which 
is homeomorphic to x 

The Reissner-Nordstrom metric g = gn e in the coordinates (t,r) is given 

by 

g^-Dde + ^dr^ + r^g^2, (85) 

where 

D^D{r) = l- — + -^ (86) 
r 

and gg2 is the standard metric on Note that the Maxwell potential A in 
these coordinates is given by 

~— dt- B cos 9 d(j) 
r 

where e = yjQ'^ -f 52, {9,(t>) G and Q, B are the electric and magnetic 
charge, respectively. 

Clearly, S0(3) acts by isometry on these spacetimes. We will refer to the 
SO(3)-orbits as (symmetry) spheres. The coordinate r is defined intrinsically 
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such that the area of the spheres of symmetry is 47rr^ (and thus should be 
thought of as a purely geometric function of the spacetime). 

One could now pose the following question: On what manifold is the metric 



(85) most naturally defined? In the above coordinates, it is clear that the metric 
component g^r is singular at r = 0, r_, r_|_ where r_, r_|_ are the roots of D. The 
computation of the curvature shows that as r — !■ the curvature blows up and 
so the singularity of r = in (851 is essential (for a very detailed description of 



these phenomena in Schwarzschild case see [22] and [29]). However, the points 
where r = , r+ form coordinate singularitietj^ which can be eliminated by 
introducing the so-called tortoise coordinate r* 

dr* (r) 1 
dr ~ D' 



We can easily see that if r_|- > r_ then 



r* {r) — r + In 



2k. 



In 



C 



where k+, k_ are given by 



K± 



r± — 7'zp 1 dD (r) 



^ri 2 dr 

In the extreme case we have r+ = r_ = M and thus k+ = k_ = 0. Then 

M2 



'{r) = r + 2Mln(r - M) - 



r- M 



C. 



(87) 



The fact that in extreme case r* is inverse linear (instead of logarithmic in the 
non-extreme case) is crucial. The constant C is taken such that r* {Q) — 0, 
where 

(88) 

The physical interpretation of the radius Q and the constants k+,k- become 




apparent in Sections 2.2 and 2.3 By introducing the coordinate system {t,r*) 
the metric becomes 

g = -Ddt^ + D {dr* f + r^g^2. (89) 



^It is the function t that is singular at these points. 
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Figure 1: The coordinate systems {t,r) and (t,r*). Each point in this diagram 
represents a sphere of symmetry. 



The tortoise coordinate r* maps the hypersurface r = r_|_ to r* = — oo. 
Therefore, although the coordinates (t, r*) capture the geometry of the region 
r > r+ in a more appropriate way than the coordinates (t, r) they still fail to 
reveal the geometric features of the neighbourhoods of the hypersurface r = r+. 
This is done by the ingoing Eddington-Finkelstein coordinates {v,r) where 

V ^t + r*. 



In these coordinates the metric is given by 

g = -Ddv^ + 2dvdr + r^g^2 . (90) 

The coordinate vector field dy is KilHng and causal in the region 

{0 < r < r_} U {r > r+} . 

In particular, dy is everywhere timelike in this region except on the hypersurfaces 

H+ = {r = r+},-H+ = {r = r_} 

where it is null. Note that in extreme Reissner-Nordstrom we have r_ = r-|_ and 
so = The region between H"*" and disappears and the Killing vector 



field dy is everywhere causal. In view of (90) and the fact that dy is tangent 
and null on 7^+ we have that the vector dy is normal to H~^. Recall that if the 
normal of a null hypersurface is Killing then the hypersurface is called Killing 
horizon. The importance of such null hypersurfaces will become apparent later. 
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Figure 2: The coordinate system {v,r) 



The radial curves v — c, where c is a constant, are the ingoing radial null 
geodesies. This means that the null coordinate vector field dr differentiates with 
respect to r on these null hypersurfaces. This geometric property of dr makes 
this vector field very useful for understanding the behaviour of waves close to 
H+. 

We see that the coordinates {v,r) extend the domain that the coordinates 
(t,r) cover. Therefore, using the coordinates let us define 

7W = (-oo,+oo) X (0,+(X)) X (91) 

Figure [2] describes the structure of M . 

Another coordinate system that partially covers Ai is the null system {u, v) 
where 

u = t — r , 
V = t + r* 



and with respect to which the metric is 

g = -Ddudv + r'^g^2 . (92) 

The hypersurfaces v = c and u = c are null and thus this system is useful for 
applying the method of characteristics or understanding null infinity. 
At last, we mention another coordinate 

t* = t + r* - r ^v~r (93) 
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Figure 3: The coordinate system {t* ,r) 

As far as the Cauchy problem for the wave equation is concerned, a possible 
initial hypersurface is given in these coordinates by {w = r} = {t* = 0}. 

A. 2 Useful Reissner-Nordstrom Computations 

If = f^dy + frdr then the current defined in Section [s] is given by 



where the coefficients in (w,r) coordinates are given by 



[drfr) _ fr_ 
2 r 



frP' 

2 ' 



2/r 1 
F^r = D (drfv) ^, Ff {drfr) 



and in {t,r*) coordinates by 



2D r 



F * , — 

, J. r^r* 



/' / 



— — , F^ — — — — 



/' fD' 



2D r' 



where /' 



_dl_ 

dr* 



A. 2.1 The Wave Operator 

The wave operator in {v,r) coordinates is 



□ „V = Ddrdrip + 2dydri^ H — dy'lp + Rdrip + i\tp, 

r 
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"sj^dr = ( ) dr, 



where R = D' + ^ and D' = Moreover, note that 

^vdv = ^ + (~~2~~~) 
2 

Vr^r = 

and observe that in degenerate black holes the right hand side of all of them 
vanishes on the horizon. Also 

Div dy = 0, Div dr = -. 

r 

The wave operator in {t,r*) coordinates is 

□gV = ^ + r-^dr, (r^a,. V)] + f V 

and 

Div 5^. = D' + — = R. 

r 

In {u,v) coordinates 

4 D' 

Ogip = -—dudy{rtp) tp + /^tp 

Ur r 

and ^ 

Div dy = — Div du = 1 • 

2 r 

Note that if a G K then 



A. 2. 2 The Non-Negativity of the Energy-Momentum Tensor T 

It is essential to know how the energy momentum tensor T depends on the 
derivatives of tp. We use the coordinate system {v, r, 9, <p) and suppose that 

Vi = V = {V\V\0,0) 

and 

^2 =n= (n^n^ 0,0). (94) 

The reason for the above notation is that everytime we apply the vector field 
F as a multiplier we have to contract with the normal n to the boundary . 
We proceed by computing 
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when V and n are causal future directed vectors. First note that under these 
assumptions we have that V" and are non-negative. Indeed, we have 

g (n, n) = -D {n^'f + 2n'"n' < 0, 
g{n,T) = -Dn'' +rf < 0. 

If n'" < then we would have 

^Dn" + 2rf >0 > -Dn" + ^ 
Dn" 

—— <rf < Dn\ 



These two inequalities imply that > which is contradictiorp^ Similarly, 
we take > 0. In view of equation (131 we have 



¥"11"+ 



2 2 2 



2 2 2 



First observe that 



\DY"n}' - iy^n'- - \vn'' = ^ (-g(n,n)) + {~g{y,Y)) 
2 2 2 4n" ^ ' 4]/« v i'v ' ji 



1 



5(^,n). 



Clearly the above makes sense even if ?i„ = since n^j g(n^n). Furthermore, 



1 



1 



1 



1 



-D^y^n'' - -DY"n'' - -DV'n'' + V^rf = -D^Vn" + 



rvv , 9iy,V)gin,n) 



2V" 2n^ 



Therefore, if we denote wy = 2{v^y^ (^.^(^i^)) ^-nd similarly for n then we 
obtain 



+ • UJn 



""Clearly if ra is timelike then > 0. 
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Note also that 



D 
2" 



D2 



+ 2u!y ■ UJn 



where the first fraction on the right hand side is well defined (even on T-L^) since 
LOy ■ can vanish at most like D^. Indeed, if > therp^ my can vanish 
on 'H'^ at most like D and similarly for n. If, in addition, n is timelike (as will 
be the case in this paper) then the denominator can vanish at most like D and 
thus the fraction vanishes on . Clearly away from the denominator is 
strictly positive. Then, 



Vn" 1 




+ 2LUy ■ UJn 



(a.^f +[w('^)](a.^f + 



D2 



+ 2LUy ■ UJn 



+ 2ujy ■ L0„ 



(95) 



Note that if V is null and tangent to then the coefficient of the transversal 
derivative vanishes. Therefore, (95) will allow us to understand the rate of the 



degeneration of this coefficient for several causal multipliers V . 



A. 3 Penrose Diagrams 

The Penrose diagram of a fundamental domain of subextreme Reissner-Nordstrom 
in the range < e < M is 



^■0 ,-' 




1+ 



I8j£ yv vanishes on then we can argue similarly bringing back the factor which is 
outside the brackets. However, in this paper we only need to consider the case where V" > 0. 
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One can glue all these solutions together to obtain the following "maximar 
solution in the rangcF^ < e < M 



The geometry of the above diagram was discovered by Graves and Brill in 
1960. Before presenting the Penrose diagram for extreme Reissner-Nordstrom, 
it is useful to consider the following subset of the above diagram 




i 



In the extreme case the regions between the two horizons collapses to a single 
hypersurface = Ti.^ . One obtains in fact: 

^^Note that in the schwarzschild case e = 0, the "maximal" solution is rather different and 
one of its features is that it is globally hyperbolic. 
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''a'._ 
» 

The maximally extended solution is 



9,/ 




o' : 



B Stokes' Theorem on Lorentzian Manifolds 

If 7?, is a pseudo-Riemannian manifold and P is a vector field on it then we have 
the identity 

Jtz Jdn 

which is an application of Stokes' theorem. Both integrals are taken with respect 
to the induced volume form. Note that ng-fi is the unit normal to dTZ and its 
direction depends on the convention of the signature of the metric. For example, 
if 7?, is a Riemannian manifold then ni)n is the outward directed unit normal to 
dTZ. On the other hand, for Lorentzian metrics with signature (— , +, +, +) the 
vector ngu is the inward directed unit normal to dTZ in case dTZ is spacelike and 
the outward directed unit normal in case dTZ is timclike. If dTZ (or a piece of 
it) is null then we take a past (future) directed null normal to dTZ if it is future 
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(past) boundary. Recall that a piece of the boundary is future (past) if the past 
(future) of points on it lie in TZ. The following diagram is embedded in M^"*"^ 



If u is a function on TZ and V a vector field and if P = uV then 
V^P^ = (uV) = (V^u) + uV^ V = Vit • V + u Div (V) , 
therefore, 

/ Vu • V + / w Div (V) = / uV ■ nan- 
If we assume that V = Vu then 

/ Vm • Vu + / u{\I\v)— / uWv ■ noTi. 
Jtz Jtz Jan 

Note that a special case of the above equality is to consider TZ = and V = y/v 
and take 

/ fu-fv = - u{^v). (96) 

Moreover, if we set V = Wgfr then 

f d \ d f d \ dv ( d 

\ ox^ J ax^ \ ox^ J ax^ \ ax^ 

yields 

n{^)^h (£) + L (^) ^ L ™^ ■ ^''^ 

Note also that another form of Stokes' theorem is 

(9^ J)dxi A dx^ A ... A dx" = f /dx^ A ... A dx". 

Indeed, we just have to consider the form w — fdx^ A ... A dx" and notice that 
du! — {dxif )dx^ A dx^ A ... A dx". Clearly, the only components of the boundary 
that contribute on the right hand side are these where {xi = c}, c constant. Of 
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course, this "less geometric" identity is completely consistent with (97). Note 
that the appearance of the divergence in (971 is due to the fact that when a 
coordinate vector field acts on the factor of the volume form then we take the 
divergence of the field times this factor, i.e. 



diy/g ^ (Diva,) 



C Sobolev Spaces 

Given a function ip : A4 ^ R, where g) is Riemannian manifold we consider 
the tensors V'0,V^V': etc. where V is the associated Levi-Civita connection. For 
example, the Hessian is given by 

Given two tensor fields Tiji^-.-ifc and Sj^j^,,,j^, of the same type, we define the 
induced inner product on tensors by 

IT '? \ — (7^1 • n*2J2 . . „ik3k .T -.<?-- 
\±p,Opi —y y ■■■ y -1-1112. ..Ik '^'Ji'j2...3ki 

where p & and similarly we define an inner product on tensor fields by 



JM 



Therefore, if T = V'^V and S = V^ct), fc e N, then we define the homogeneous 
Sobolev inner product by 

JM 

Note that this is an inner product provided t^, (j) either decay at infinity of vanish 
at the boundary of M. . The Sobolev inner product is defined by 



and thus the Sobolev space H'^{M) is the set of functions such that the Sobolev 
norm 



is finite. Note that if ip is not sufficiently regular then the derivatives in 



are defined in a weak (distributive) way using (97). Then the Sobolev space 
H^{M) is complete and thus Hilbert space. In other words, -0 € H^{M) if and 
only if all the derivatives up to the k-th order are in L?'{M). Finally, we define 
the space H^^^{M.) to be the set of all functions such that the above expression 
is finite over all compact subsets of M . 
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Now according to the Sobolev inequality, if n is the dimension of the space 
and q> 1 and k < ^ such that 



1 



1 



k 



2 



n 



then 



Moreover, if A; > | , then 



IIV'IIl^. <C||^| 



where the Sobolev constants depend only on the geometry of the space and on 
k,q,n. 

D Elliptic Estimates on Lorentzian Manifolds 

Let us suppose that {A4,g) is a globally hyperbolic time-orientable Lorentzian 
manifold which admits a Killing vector field T. We also suppose that Ai is 
foliated by spacelike hypersurfaces S,-, where = (j)^ (So)- Here, Eg is a 
Cauchy hypersurface and <pr is the flow of T. 

Let A'' be a ^T-invariant timelike vector field and constants -Bi,B2 such 
that 



We will first derive the required estimate in Sq which for simplicity we 
denote by S. For each point p G E the orthogonal complement in TpAi of the 
line that contains N is 3-dimensional and contains a 2-dimensional subspace of 
the tangent space TpS. Let X2,X3 be an orthonormal basis of this subspace. 
Let now Xi be a vector tangent to S which is perpendicular to the plane that 
is spanned by X2,X3. Note that the line that passes through Xi is uniquely 
determined by N and E. Then, the metric g can be written as 



with respect to the frame {N,Xi,X2,X^). Then, if we set \g\ = g^N " dXiXi — 



< Bi < -g{N,N) < B2. 



9 



/ gNN gNXi 0\ 
gNXi 9x1X1 



10 
\ 1/ 



9nXi' the inverse is 



9 



-1 



/ \^\9XiXi -j^ffivxi 0\ 
--^\9NXi ]J\9nn 



10 

V 1/ 
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Let ft-s be the induced Riemannian metrie on the spacelike hypersurface S. 
Clearly, in general we do not have /i^ = g^-' ■ Indeed 

^9x,x, 0^ 
/lE = I 10 
ly 



and therefore, 



0^ 

/is = I 10 



ly 

Let : — >■ M satisfy the wave equation. Then, 

□gV = trg (Hess V) = ff"'^ (V^V)^^ = 

= 5°' ((vV)o, + (v^^),o)+5^^'(vV),,. 

We will prove that the operator 
is strictly elliptic. Indeed, in view of the formula 
the principal part cr of P is 
If ^ e T*E, then 

= t—\9nn£,i + CI + ^3 

= m\, 

where the ellipticity constant 6 > depends only on S. Moreover, if i/' satisfies 
□gV = then 



|2 

'/30 



112 



where C is a uniform constant that depends only on the geometry of E and the 
precise choice of N. Therefore, if tp can be shown to appropriately decay at 
infinity then by a global elliptic estimate on S we obtain 



Hi(E) 



+ 



2 



for some uniform positive constant C. 

In case our analysis is local and thus we want to confine ourselves in a 
compact submanifold E of S then by a local elliptic estimate on E we have 



Hi(E) 



< 



One can also estimate spacetime integrals by using elliptic estimates. Indeed, 
if TZ (0, r) is the spacetime region as defined before, then 



/ / |Vw| dg^ -- 
Jn{o,T) Jo 



fdg^^ ] dt, 



where the integrals are with respect to the induced volume form and u : M. ^M. 

is such that u{p) = r iff j> € St-. Then Vu is proportional to and since 
T{u) = 1, Vu is (/),--invariant. Therefore, |Vw| is uniformly bounded. If now / 
is quadratic on the 2-jet of tp then 



TC(0,t) 



fdgn 



df 



< C 

c< 



/•R.(0,t) 

In applications we usually use these results away from "H"*" where we commute 
with T and we use the X estimate. We can also use this estimate even if S 
(and TV) crosses T-L^ , provided we have commuted the wave equation with N 
and A^A'' (recall that we need commutation with NN only for degenerate black 
holes). 
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